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ABSTRACT 


A  relatively  simple  procedure  has  been  developed  to  include  variable  boundary 
layer  edge  entropy  in  an  inviscid  flow  solution  on  the  windward  ray  of  axisym- 
metric  blunt  bodies.  Variable  boundary  layer  edge  entropy  is  significant  down¬ 
stream  of  the  nose  region  on  blunt  bodies  due  to  entropy  gradients  along  the 
curved  bow  shock.  Laminar  and  turbulent  boundary  layer  growth  models  are  used 
in  a  mass  balancing  technique  to  account  for  the  effects  of  entropy  gradients 
produced  at  the  bow  shock.  The  solution  is  driven  by  the  inviscid  surface  pres¬ 
sure  distribution,  mass  flow  within  the  boundary  layer,  and  bow  shock  shape.  The 
code  takes  about  10  min/case  to  run  on  an  IBM  AT  personal  computer.  Heat 
transfer  comparisons  to  experimental  data  and  other  heating  rate  codes  are  favor¬ 
able.  The  method  is  accurate  from  M«,  -  2.8  to  15  for  cones  with  flow  inci¬ 
dence  angles  up  to  25°. 
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1.0  INTRODUCTION 


Predicting  aerothermal  heating  along  the  windward 
meridian  of  axisymmetric  blunt  bodies  in  supersonic  and 
hypersonic  flight  is  a  common  engineering  problem.  A 
simple,  fast,  and  inexpensive  way  to  calculate  heating  is 
to  correlate  boundary  layer  (BL)  edge  conditions  to  a 
heat  transfer  coefficient.  Inherent  to  the  method  is  deter¬ 
mining  the  flow  properties  at  the  BL  edge  through  the 
solution  of  greatly  simplified  forms  of  the  Navier-Stokes 
equations.  Ever  present  is  the  drive  toward  simpler  meth¬ 
ods  that  sacrifice  as  little  accuracy  as  possible  when  com¬ 
pared  to  complete  Navier-Stokes  solutions  to  the  body 
surface. 

The  classical  (and  simplest)  approach  for  calculating 
bL  edge  conditions  is  to  neglect  coupling  between  the 
inviscid  flow  region  and  the  BL,  thereby  treating  the 
BL  edge  as  an  isentropic  surface.  The  properties  along 
the  BL  edge  are  found  by  following  the  stagnation 
streamline  as  it  passes  through  the  normal  bow  shock 
and  then  by  assuming  an  isentropic  expansion  from  the 
stagnation  point  around  the  body  surface.  The  method 
requires  only  the  pressure  distribution  and  normal  shock 
entropy  for  a  complete  thermodynamic  definition  of  the 
BL  edge.  Currently,  a  computer  code  named  ABBOHT 
(Axisymmetric  Blunt  Body  Heat  Transfer)  is  used  to  cal¬ 
culate  cold  (isothermal)  wall  heat  transfer  to  the  wind¬ 
ward  ray  of  blunted  bodies  in  just  such  a  manner. 

In  1953  Ferri  and  Libby1,2  found  that  entropy  gra¬ 
dients  induced  by  bow  shock  curvature  around  a  blunted 
body  produce  a  coupling  effect  between  the  inviscid  re¬ 
gion  and  the  BL  that  can  invalidate  the  classical  solu¬ 


1  A.  Ferri  and  P  A.  Libby,  “Note  on  the  Interaction  Between  the 
Boundary  Layer  and  the  Inviscid  Flow,”  J.  Aeronaut.  Set.  21, 
130  (1954) 

'A.  Ferri,  “Some  Heat  Transfer  Problems  in  Hypersonic  Flow,” 
in  Aeronautics  and  Astronautics,  Pergamon  Press,  New  York,  pp. 
344-377  (i960) 


tion.  First,  the  entropy  along  the  BL  edge  will  continu¬ 
ously  vary  aft  of  the  nose  region;  second,  the  velocity 
gradient  at  the  predicted  BL  edge  ( du/dy )  will  not  be 
zero  (producing  nonsimilar  BL  growth)  as  previously 
assumed.  The  latter  effect  is  significant  when  the  ex¬ 
ternal  stream  vorticity  is  of  the  same  order  of  magnitude 
as  the  average  BL  vorticity  (i.e.,  low  Reynolds  no.  com¬ 
bined  with  high  Mach  no.  flow).  External  flow  vorticity 
effects  are  small  at  high  Reynolds  no.  flow  and  on  slen¬ 
der  blunted  cones  (half  angles  less  than  25°).  Since  that 
is  the  flow  regime  commonly  used  in  ABBOHT,  it  was 
desired  to  upgrade  the  code  to  include  variable  entro¬ 
py  effects  while  accepting  the  BL  similarity  limitations. 

Modifications  to  the  classical,  first-order  approach  to 
include  shock  curvature  effects  were  investigated  from 
1960  to  1968  by  Zakkay  and  Krause,  Wilson,  Rotta  and 
Zakkay,  and  Murzinov.36  All  found  the  need  to  ac¬ 
count  for  variable  BL  entropy  through  “entropy  layer 
swallowing”  or  “streamline  swallowing.”  That  event 
occurs  as  follows:  BL  growth  along  a  body  is  caused 
by  mass  entrainment  of  the  inviscid  external  flow  with 
streamlines  that  have  passed  through  an  increasingly 
oblique  portion  of  the  bow  shock.  Mass  entrainment 
eventually  swallows  the  high-entropy  BL  from  the  nose 
region.  Consequently,  the  BL  downstream  of  the  nose 
region  consists  of  lower  entropy  streamlines  resulting 
in  higher  BL  velocities  and  higher  heat  transfer  rates 
than  for  the  classical  method.  Low  BL  velocities  and 
heat  transfer  at  aft  body  stations  in  the  old  version  of 
ABBOHT  are  documented  in  Ref.  7. 


3V.  Zakkay  and  E.  Krause,  “Boundary  Conditions  at  the  Outer 
Edge  of  the  Boundary  Layer  on  Blunted  Conical  Bodies,”  AIAA 
J.  1,  1671-1672  (1963). 

4R.  E.  Wilson,  “Laminar  Boundary  Layer  Growth  at  Hypersonic 
Speeds,”  J.  Spacecr.  Rockets  2,  490-496  (1965). 

SN.  R.  Rotta  and  V.  Zakkay,  “Effects  of  Nose  Bluntness  on  the 
Boundary  Layer  Characteristics  of  Conical  Bodies  at  Hypersonic 
Speeds,”  Astronaut.  Acta  13,  507-516  (1968). 

61.  N.  Murzinov,  “Laminar  Boundary  Layer  on  Blunt  Bodies,  Al¬ 
lowing  for  Vorticity  of  the  External  Stream,”  NASA  TT  F-l  1007 
(Jun  1966). 

7J.  B.  Kouroupis,  “Comparison  of  Aerodynamic  Heating  Between 
Sharp  and  Blunt  Bodies,”  JHU/APL  BBE/EM-5292  (12  Mar 
1986) 
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2.0  STREAMLINE  SWALLOWING  METHODS 


2.1  MASS  BALANCING 

Two  simplified  (first-order)  approaches  have  been  de¬ 
veloped  to  account  for  entropy  layer  swallowing  in  an  in- 
viscid  code.  The  first  is  called  mass  balancing,  whereby 
the  mass  flow  in  the  BL  at  the  station  of  interest  is 
matched  with  the  mass  flow  in  a  free-stream  cylinder  of 
appropriate  radius  (ysh )  as  shown  in  Fig.  1,  which  is 
based  on  a  figure  in  Ref.  8.  The  value  of  ysh  is  defined 
as 

m 

ysh  ~  /  • 

V/TTpooW*. 


Figure  1  Streamline  swallowing  on  a  blunted  body. 

The  bow  shock  angle  at  ysh  is  then  used  to  calculate 
the  entropy  behind  the  shock,  which  is  used  with  the 
inviscid  surface  pressure  to  determine  BL  edge  condi¬ 
tions..  Thus,  the  bow  shock  shape,  inviscid  pressure  dis¬ 
tribution,  and  BL  mass  flow  are  required  for  the  solu¬ 
tion  of  the  above  equation. 

In  1970-71  Mayne,  Adams,  and  Dyer8,9  built  on  the 
work  done  by  Patankar  and  Spalding8 9 10  to  develop 


8A.  W.  Mayne,  Jr.  and  J.  C.  Adams,  Jr.,  “Streamline  Swallow¬ 
ing  by  Laminar  Boundary  Layers  in  Hypersonic  Flow,”  AEDC- 
TR-71-32,  Arnold  Engineering  Development  Center,  Tullahoma, 
Tenn.  (Mar  1971). 

9A.  W.  Mayne,  Jr.  and  D.  F.  Dyer,  “Comparisons  of  Theory  and 
Experiment  for  Turbulent  Boundary  Layers  on  Simple  Shapes  at 
Hypersonic  Conditions,”  Proc.  1970  Heat  Transfer  and  Fluid 
Mechanics  Institute,  Stanford  University  Press,  168-188  (1970). 

I0S.  V.  Patankar  and  D.  B.  Spalding,  “A  Finite-Difference  Proce¬ 

dure  for  Solving  the  Equations  of  the  Two-Dimensional  Bound¬ 

ary  Layer,”  Ini  J.  Heat  Mass  Transfer  10,  1389-1411  (1967). 


heat  transfer  codes  tint  use  mass  balancing  to  account 
for  streamline  swallowing  for  both  laminar  and  turbu¬ 
lent  BL  flow.  Here,  BL  equations  are  expressed  in  a 
normalized  Von  Mises  coordinate  system  and  are  solved 
using  an  implicit,  finite-difference  marching  technique. 
BL  mass  flow  is  readily  derived  from  that  solution.  A 
Method  of  Characteristics  solution  is  used  for  the  in¬ 
viscid  region  to  calculate  bow  shock  shape  and  surface 
pressure.  Mayne’s  laminar  method  accounts  for  non¬ 
similar  BL  growth,  but  the  turbulent  method  does  not. 

In  1972  Price  and  Harris11  used  a  technique  similar 
to  Mayne’s  but  solved  the  nonsimilar  BL  equations  for 
laminar,  transitional,  and  turbulent  flow.  Thus,  all  three 
BL  flow  regimes  can  account  for  the  effects  of  vorticity 
interaction.  As  with  Mayne’s  method,  an  inviscid  flow 
solution  from  the  bow  shock  to  the  body  surface  was 
added  to  define  the  surface  pressure  distribution  and 
bow  shock  shape,  thereby  allowing  the  inclusion  of  mass 
balancing.  The  program  was  revised  in  1982 12  to  in¬ 
crease  its  accuracy  and  computational  efficiency. 

In  1974  DeJarnette  and  Hamilton13  devised  a  differ¬ 
ent  approach  for  using  mass  balancing  to  account  for 
streamline  swallowing.  First,  a  modified  Newtonian  sur¬ 
face  pressure  distribution  is  calculated  over  the  body. 
Inviscid  surface  streamlir.  s  are  then  defined  using  Eu¬ 
ler’s  equation  along  with  two  stream  functions.  This  is 
an  inverse  method,  where  the  bow  shock  shape  can  be 
found  from  stream  function  values.  Those  values  also 
define  the  BL  mass  flow  rate,  which  is  matched  to  free- 
stream  flow  rates.  The  method  is  iterative,  does  not  ac¬ 
count  for  nonsimilar  BL  growth,  and  is  generally  less 
accurate  than  Mayne’s  method;  it  is,  however,  simpler 
and  can  be  extended  easily  to  three-dimensional  geom¬ 
etries.  Modifications  to  make  the  basic  technique 
noniterative  were  made  by  Fivel. 14  Additional  changes 
to  make  it  simpler  (albeit  less  accurate)  for  the  rapid 


"S.  M.  Price  and  J.  E.  Harris,  “Computer  Program  for  Solving 
Compressible  Nonsimilar  Boundary-Layer  Equations  for  Lami¬ 
nar,  Transitional,  or  Turbulent  Flows  of  a  Perfect  Gas,”  NASA 
TM  X-2458  (1972). 

I2J.  E.  Harris  and  D.  K.  Blanchard,  “Computer  Program  for  Solv¬ 
ing  Laminar,  Transitional,  or  Turbulent  Compressible  Boundary- 
Layer  Equations  for  Two-Dimensional  and  Axisymmetric  Flow,” 
NASA  TM-83207  (Feb  1982). 

"F.  R.  DeJarnette  and  H.  H.  Hamilton,  “Aerodynamic  Heating 
on  3-D  Bodies  l  tcluding  the  Effects  of  Entropy  Layer  Swallow¬ 
ing,”  J.  Spaced \  Rockets  12,  5-12  (Jan  1975). 

IJH.  J.  Fivel,  “Numerical  Flow  Field  Program  for  Aerodynamic 
Heating  Analysis,”  AFFDL-TR-79-3128,  Vol.  I,  Air  Force  Flight 
Dynamics  Laboratory,  Wright  Patterson  Air  Force  Base,  Ohio 
(1979). 
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and  economical  computation  of  heating  rates  and  sur¬ 
face  temperatures  throughout  an  entire  trajectory  were 
made  by  DeJarnette  et  al.15 

2.2  INVISCID  APPROACH 

The  second  simplified  approach  used  to  calculate  vari¬ 
able  BL  edge  entropy  was  born  of  the  need  to  account 
for  entropy  layer  swallowing  over  three-dimensional 
shapes  (i.e.,  space  shuttles),  where  mass  balancing  tech¬ 
niques  become  very  complicated,  and  at  high-incident 
angles,  where  mass  balancing  using  similar  BL  develop¬ 
ment  overpredicts  heat  transfer.  This  approach  accounts 
for  entropy  layer  swallowing  by  assuming  BL  edge  con¬ 
ditions  to  be  equal  to  inviscid  properties  at  the  hypothet¬ 
ical  BL  thickness.  The  BL  thickness  is  found  through 
its  relation  to  the  momentum  thickness,  which  is  de¬ 
rived  from  the  momentum  integral  equation  using  an 
expression  relating  Reynolds  no.  and  skin  friction.  A 
complete  inviscid  flow  solution  is  required  for  this 
method. 

The  inviscid  technique  was  originally  developed  by 
Zoby,  Moss,  and  Sutton16  in  1980.  BL  edge  conditions 
found  in  this  manner,  when  combined  with  the  approx¬ 
imate  heating  relations  given  in  Ref.  16,  yield  favorable 
heat  transfer  comparisons  over  blunted  bodies  for  a  va¬ 
riety  of  wind  tunnel  and  flight  conditions,  although  the 
technique  has  been  used  primarily  for  high-incident- 
angle  and  three-dimensional  problems. 


2.3  NEW  ABBOHT  METHOD 

The  method  chosen  for  incorporation  into  ABBOHT 
is  a  simplified  mass  balancing  version.  It  retains  the  ad¬ 
vantage  of  requiring  only  the  invite!  J  pressure  distri¬ 
bution  'as  with  DeJamette’s  method)  instead  of  the 
entire  inviscid  flow  field  (as  with  the  Mayne,  Price  J 
Zoby  methods).  Moreover,  it  is  simpler  than  either  of 
the  mass  balancing  methods  discussed  above  because  it 
neither  solves  the  complex  BL  equations  (as  with  the 
Mayne  and  Price  methods)  nor  requires  the  use  of  mul¬ 
tiple  stream  functions  (as  with  DeJarnette’s  method). 
This  new  method  uses  a  unique  “tangent-blunted-cone” 
pressure  technique  (already  contained  in  the  old  ver¬ 
sion  of  ABBOHT)  conceived  at  APL  by  L.  B.  Weck- 
esser  and  developed  by  L.  L.  Perini. 17  This  approach 
does  not  include  nonsimilar  BL  growth,  so  its  ability  to 
do  large  half-angle  cones  is  limited.  For  computational 
efficiency,  only  simple  relations  are  used  for  BL  growth 
models  and  heating  rate  correlations. 

Presented  below  is  a  description  of  the  new  method 
to  calculate  streamline  swallowing  effects  as  incorpo¬ 
rated  into  ABBOHT.  The  new  method  also  is  compared 
to  the  old  ABBOHT  version;  experimental  data;  the  pre¬ 
diction  methods  of  Mayne,  DeJarnette,  and  Zoby;  and 
viscous  shock  layer  (VSL)  predictions. 


3.0  DISCUSSION 


The  first  step  in  calculating  BL  edge  conditions  in 
ABBOHT  is  to  determine  the  inviscid  surface  pressure 
distribution  as  described  in  Ref.  17.  To  use  the  mass 
balancing  technique  to  account  for  streamline  swallow- 

i5F.  R.  DeJarnette,  L.  A.  Kania,  and  A.  Chitty,  “Aerodynamic 
Heating  and  Surface  Temperatures  on  Vehicles  for  Computer- 
Aided  Design  Studies,”  A1AA  Paper  83-0411  (Jan  1983). 

I6E  V.  Zoby,  J.  N.  Moss,  and  K.  Sutton,  “Approximate  Convec¬ 
tive  Heating  Equations  for  Hypersonic  Flows,”  J.  Soacecr.  Rockets 
18,  64-70  (Jan  1981). 


ing,  the  inviscid  pressure  distribution  is  coupled  with 
a  BL  (mass  flow  rate).  To  incorporate  the  required  pro¬ 
cedure  into  ABBOHT,  two  geometrical  additions  must 
first  be  made:  laminar  and  turbulent  BL  growth  models 
for  compressible  flow  (ABBOHT  does  not  include  tran¬ 
sition  region  heating);  and  the  bow  shock  shape.  A  dis- 


,7L.  L.  Perini  and  J.  B.  Kouroupis,  “A  Simplified  Procedure  for 
Calculating  Pressure  Distnbutions  Over  Axisymmetnc  Blunt  Nosed 
Bodies,”  JHU/APL  BBE/EM-5230  (14  May  1985). 
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cucsion  of  how  the  additional  ABBOHT  constituents 
are  calculated  is  presented  first,  followed  by  an  over¬ 
view  of  the  complete  procedure  for  calculating  heat 
transfer  rates  and,  finally,  a  comparison  of  results.  The 
ABBOHT  code  is  written  in  Fortran  77  and  was  run 
on  an  IBM  AT  desktop  computer  for  all  results. 


3.1  BOUNDARY  LAYER  GROWTH 

3.1.1  Laminar  Boundary  Layer 

The  governing  equation  used  for  mass  flow  rate  in 
a  BL  over  a  conical  body  is 

tit  =  2itrpeue(b-b*)  . 


The  resulting  equation  for  mass  flow  rate  is 


2irrue peSTe 
V3  v JueS/vw  Tw 


[2.007  p  +  2.01  [l  +  ^  Mg 
-  2.003  A/f2]  . 


To  account  for  the  initial  BL  thickness  caused  by  the 
blunted  nose,  the  stagnation  BL  thickness  is  first  cal¬ 
culated  from  Ref.  19  using 


Therefore,  5  and  b*  must  be  modeled.  Since  it  is  desired 
to  know  the  momentum  thickness  Reynolds  no.  as  a 
criterion  for  BL  transition,  6  also  must  be  calculated. 
All  BL  growth  calculations  start  at  the  nose/body  tan- 
gency  point.  Ness18  modeled  the  three  BL  parameters 
for  compressible,  laminar  flow  over  a  sharp  cone  by 
solving  the  energy  and  X-momentum  equations  assum¬ 
ing  Pr  -  1,  Pe,  Te,  Me,  ue,  Tw,  vw>  7  =  constant 
and  0  =  0.  While  these  are  crude  assumptions  for  a 
blunted  cone  (especially  in  the  expansion-compression 
region  at  the  nose/body  tangency  point),  they  have  been 
demonstrated  herein  to  be  adequate  for  use  in  this  mod¬ 
el.  The  resulting  equations  for  5,  6*,  and  6  are 


The  equation  (obtained  from  Ref.  20)  for  the  stagna¬ 
tion  point  velocity  gradi^m  is  empirical  and  valid  for 
all  free-stream  Mach  numbers. 

The  value  of  bs  is  then  added  to  the  sharp  cone  BL 
thickness  at  the  nose/ body  tangency  point,  yielding  the 
initial  BL  thickness  at  the  start  of  the  calculations. 


S  Te 

V3  \keS/  vw  Tw 

3.00  —  +  2.01  [l  +  —  Me2 
Te  L  2 

_  i.62  nl  A/f2]  , 

S  Te 

V3  'JueS/vw  Tw 


3.1.2  Turbulent  Boundary  Layer 

Whereas  the  laminar  BL  equations  lend  themselves 
to  a  closed-form  solution  for  BL  growth,  the  turbulent 
equations  do  not.  Many  semi-empirical  and  correlation 
methods  to  predict  turbulent  BL  development  are  avail¬ 
able.  The  method  chosen  for  use  in  ABBOHT  (the  new 
version)  is  a  flat  plate  formu'ation  originally  proposed 
for  incompressible,  turbulent  flow  by  Head21  in  1958 
and  then  extended  for  compressible  flow  by  Green22 
in  1967.  This  method  was  chosen  for  its  proven  accuracy 
in  the  flow  regime  of  interest  and  because  it  solves  direct¬ 
ly  for  the  quantity  (b  -  §*),  which  is  used  for  model¬ 
ing  the  BL  mass  flow  rate. 


0.993  —  +  0.383  —  A/,2!  , 
Te  2  ) 


0.383  S  Te 
V3  \fu~S/vH  Tw 


'V.  Ness,  “Dynamics  of  Viscous  Fluids,”  Course  Notes,  Dept,  of 
Aerospace  Engineering,  West  Virginia  University  (Sep  1965) 


,9R.  W.  Truitt,  Fundamentals  of  Aerodynamic  Heating,  The  Ronald 
Press  Co.,  New  York  (I960). 

^L.  L.  Perini,  "Compilation  of  Experimental  Stagnation  Point  Ve¬ 
locity  Gradients  and  Heat  Transfer  Data  in  Subsonic  and  Super- 

i  sonic  Flow,”  JHU/APL  AEO-75-29  (14  Aug  1975). 

21M.  R.  Head,  “Entrainment  in  the  Turbulent  Boundary  Layer,” 
Aeronautical  Research  Council  Reports  and  Memoranda,  R  &  M 
No.  3152  (Sep  1958). 

22J.  E.  Green,  “The  Prediction  of  Turbulent  Boundary  ^ayer  De¬ 
velopment  in  Compressible  Flow,”  J.  Fluid  Mech.  31,  753-778 
(1968). 
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The  underlying  assumption  for  Head’s  method  is  that 
the  BL  entrainment  rate  is  proportional  to  the  velocity 
defect  in  the  outer  portion  of  the  BL.  The  velocity  de¬ 
fect  can  be  approximated  as  the  function  of  a  modi¬ 
fied  form  parameter  (//, )  multiplied  by  the  external 
stream  velocity,  or 

d(ue  A) 

entrainment  rate  -  -  =  ueF(Hx)  , 

dx 


where 

Hx 

A 

and 

A 

II 

o> 

1 

09 

# 

It  was  empirically  determined  that 

F  =  0.0306  (Hx  -  3.0)  ~0  653 
for  incompressible  flow. 

Using  Head’s  assertion  that  density  variation  in  the 
outer  part  of  the  BL  had  little  effect  on  the  entrain¬ 
ment  mechanism,  Green  inferred  that  the  incompress¬ 
ible  relationship  between  F  and  Hx  also  is  valid  for 
compressible  flow.  However,  Head’s  relationship  be¬ 
tween  H  (which  is  used  in  the  momentum  integral  equa¬ 
tion  below)  and  Hx  was  inaccurate  for  compressible 
flow;  therefore,  Green  provided  empirical  transforma¬ 
tion  functions  between  the  compressible  and  incompress¬ 
ible  versions  of  the  form  parameters. 

Since  F  is  assumed  to  be  unaffected  by  density  vari¬ 
ation,  Green  accounted  for  compressibility  effects  simply 
by  inserting  density  into  the  entrainment  rate  equation, 
yielding 

d 

—  (p<.MeA)  =  peueF  . 

dx 

By  expanding  the  derivative  and  substituting 


M\  due  1  dpe 

- —  for  —  —  , 

ue  dx  pe  dx 

Green  obtained  the  following  compressible  entrainment 
rate  equation: 


dA 

dx 


,  A  due 

F  +  (Ml  -  1)  -  ~  * 
u„  dx 


The  momentum  integral  equation  for  compressible  flow 
over  a  flat  plate  is 


dd  cf 

dx  2 


(H  +  2  -  Ml) 


e 

ue 


due 

dx 


The  conventional  form  parameter,  //(for  compress¬ 
ible  flow),  in  the  above  equation  is  found  through  an 
empirical  transformation  function  together  with  an  as¬ 
sumed  quadratic  BL  temperature  distribution,  which  is 
good  for  isothermal  walls  and  small  pressure  gradients. 
The  compressible  skin  friction  (cy)  is  dcu  rmined  using 
Eckert’s  reference  temperature  transformation  between 
incompressible  and  compressible  flow.  Since  this  growth 
model  was  derived  for  flat  plates,  the  BL  growth  rate 
was  reduced  by  a  factor  of  2.25° 2  to  compensate  for 
axisymmetric  effects.  The  correction  factor  for  turbu¬ 
lent  flow  was  derived  using  the  Mangier  transforma¬ 
tion  (see  Appendix  C  of  Ref.  23).  Three-dimensional 
effects  on  BL  growth  are  not  accounted  for,  effective¬ 
ly  yielding  a  sharp  core  BL  growth  model.  As  with  the 
laminar  BL  model,  the  sharp  cone  turbulent  growth 
model  has  been  demonstrated  herein  to  be  adequate  for 
use  in  this  heat  transfer  method.  It  is  not  the  intent  of 
the  work  reported  here  to  accurately  model  the  struc¬ 
ture  of  the  BL,  but  only  to  estimate  its  mass  flow  rate. 

By  marching  the  entrainment  and  momentum  integral 
equations  using  an  explicit,  first-order-accurate  tech¬ 
nique,  the  development  of  the  necessary  parameters  (0 
and  A)  for  mass  flow  rate  predictions  and  the  momen¬ 
tum  thickness  Reynolds  no.  is  readily  accomplished. 
Turbulent  BL  mass  flow  rate  is  simply 

m  =  2zrpeueA  . 


3.2  BOW  SHOCK  SHAPE 


The  bow  shock  shape  equation,  empirically  derived 
by  Billig,24  is  used  for  this  model.  Billig  assumes  a 
hyperbolic  shock,  asymptotic  to  the  sharp  cone-attached 
shock  angle.  Since  only  the  bow  shock  angle  is  required 
for  streamline  swallowing,  the  equation  for  coordinates 
of  the  shock  given  by  Billig  was  (’ifferentiated  to  yield 
the  slope  as 


i\r 

tan  X  =  — 

ySh 


y\h  tan2??^05 
>  * 


where 

/  0.54  \ 

Rc  =  1.143  Rn  exp  ( - ^  ) 

c  "  \(MX  -  l)'2/ 


and  r)  is  the  sharp  cone-attached  shock  angle  obtained 
from  Simon  and  Walter25  as 


23C.  Gazley,  Jr.,  “Theoretical  Evaluation  of  the  Turbulent  Skin- 
Friction  and  Heat  Transfer  on  Cone  in  Supersonic  Flight,” 
^  R49A0524,  General  Electric  Co.  (Nov  1949). 

24F.  S.  Billig,  “Shock-Wave  Shapes  around  Spherical-  and 
^Cylindrical-Nosed  Bodies,”  J.  Spacecr.  Rockets  4, 822-823  (1967). 
2!W.  E.  Simon  and  L.  A.  Walter,  “Approximations  for  Superson¬ 
ic  Flow  Over  Cones,”  A1AA  J.  1,  1696-1697  (1963). 
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sin  i]  = 


+ 


7+1 

~2~ 


Thus,  for  a  given  ysh,  Rn,  and  a,  the  slope  of  the  bow 
shock  can  be  easily  found. 


4.0  PROCEDURE 


4.1  LAMINAR  AND  TURBULENT  FLOW 

Streamline  swallowing  calculations  begin  at  the  nose/ 
body  tangency  point.  For  laminar  BL  flow,  the  proce¬ 
dure  progresses  as  foilows: 


1 .  Step  surface  distance. 

2.  Calculate  new  surface  pressure  (constant  down¬ 
stream  of  nose  region). 

3.  Use  current  value  of  BL  edge  entropy  with  new 
pressure  to  estimate  new  BL  edge  properties. 

4.  Calculate  new  mass  flow  rate. 

5.  Calculate  new  radius  of  free-stream  cylinder  to 
match  mass  flow  rate.. 

6.  Calculate  new  bow  shock  angle  at  given  ysh. 

7. ;  Calculate  new  entropy  on  downstream  side  of 

shock. 

8.  Compare  entropy  from  step  7  with  previous  value. 

9.  If  not  converged,  use  entropy  from  step  7  and  re¬ 
peat  procedure  beginning  with  step  3  until  con¬ 
verged.  All  BL  edge  properties  are  now  thermody¬ 
namically  consistent  and  the  process  repeats  again 
at  step  1  until  the  transition  Reynolds  no.  is 
reached. 


The  quantities  <5  and  6*  are  continually  calculated  for 
laminar  flow  and  their  last  two  values  are  saved  so  that 
the  turbulent  equations  may  be  pr  oerly  initialized.  Af¬ 
ter  calculating  initial  values  for  .  turbulent  entrain¬ 
ment  and  momentum  equations,  ine  turbulent  stream¬ 
line  swallowing  procedure  continues  as  follows:. 


3.  Step  entrainment  and  momentum  equations  us¬ 
ing  backward  differencing  to  find  new  values  for 
A  and  0., 

4.  Use  new  values  for  A  and  0  with  BL  edge  condi¬ 
tions  to  calculate  new  form  parameter  and  mass 
flow  rate. 

5.  Calculate  new  radius  of  free-stream  cylinder  to 
match  mass  flow  rate. 

6.  Calculate  new  bow  shock  angle  at  given  ysh. 

7.  Calculate  new  entropy  on  downstream  side  of 
shock. 

8.  Use  new  entropy  with  new  surface  pressure  from 
step  2  to  calculate  new  BL  edge  conditions. 

9.  Compare  entropy  from  step  7  with  previous  value. 

10.  If  not  converged,  use  entropy  from  step  7  and  re¬ 
peat  procedure  beginning  with  step  4  until  con¬ 
verged. 

1 1 .  Repeat  beginning  with  step  1  until  end  of  body 
is  reached. 


There  are  three  parameters  driving  the  streamline 
swallowing  solutions;  i.e.,  the  bow  shock  angle,  the  body 
radius,  and  the  surface  distance.  Consequently,  as  the 
bow  shock  angle  approaches  that  for  a  sharp  cone  well 
downstream  of  the  nose,  the  BL  edge  properties  also 
will  approach  those  for  a  sharp  cone.  This  phenome¬ 
non  is  illustrated  in  Section  5.0. 


4.2  HEAT  TRANSFER 

Many  methods  are  available  to  calculate  heat  trans¬ 
fer  rates  given  BL  edge  properties  and  wall  tempera¬ 
ture..  Therefore,  to  understand  better  the  heat  transfer 
comparisons,  a  brief  review  of  the  heating  methods  used 
by  ABBOHT  is  presented  below. 


1.  Step  surface  distance, 

2.  Calculate  new  surface  pressure. 
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4.2.1  Stagnation  Point  Heating 

Stagnation  point  heating  is  accomplished  with  Lees’s 
equation,26  modified  to  include  wall  temperature  ef¬ 
fects  by  using  Eckert’s  reference  temperature  to  evalu¬ 
ate  density  and  viscosity.  The  resulting  equation 
(assuming  Pr  =  0.72)  is 

Qs  =  0.88  ((p.m.)s  &)°5  (K  -  K)s 

4.2.2  Laminar  Heating 

Laminar  heating  also  is  accomplished  with  Lees’s 
method,  which  requires  the  stagnation  point  heating  rate 
as  an  input.  The  cold  wall  heating  rate  is  defined  as 

•  =  Qs  Hi  Ue  r 

2  V(p.  n.)s  0S  s~ 


4.2.3  Turbulent  Heating 

The  solution  of  the  momentum  integral  equation  with 
an  assumed  1/7-power  velocity  profile  is  used  for  tur¬ 
bulent  heating  relations.  Eckert’s  reference  enthalpy 
method27  is  used  to  account  for  compressibility  effects 
and  the  Reynold’s  analogy  factor  of  Colburn  is  used 
to  convert  skin  friction  into  a  heat  transfer  coefficient. 
To  convert  the  flat  plate  relation  into  one  for  an  axi- 
symmetric  body,  a  correction  factor  of  2.25° 2  was 
used,  effectively  treating  three-dimensional  heating  ef¬ 
fects  as  constant.  The  resulting  heating  rate  is 


Qw 


0.0348  p.  ue 
Re?  2  Pr 067 


(hr  -  K)  . 


5.0  COMPARISON  OF  RESULTS 


To  illustrate  the  effects  of  streamline  swallowing  on 
BL  edge  conditions,  results  from  the  old  and  new  ver¬ 
sions  of  ABBOHT  will  be  compared.  Figure  2  shows 
the  entropy  distribution  along  the  edge  of  laminar  and 
turbulent  BLs  up  to  an  S/Rn  of  300,  with  the  new 
model  compared  to  the  constant  entropy  value  used  in 
the  old  model.  Also  shown  for  comparison  is  the  sharp 
cone  BL  edge  entropy.  The  BL  edge  entropy  with  stream¬ 
line  swallowing  is  seen  to  decrease  rapidly  toward  the 
sharp  cone  asymptote  for  turbulent  BL  flow,  and  more 
gradually  for  laminar  BL  flow.  The  sharp  cone  value 
is  16%  below  the  normal  shock  value. 

Figure  3  compares  the  edge  Mach  nos.  for  the  same 
BL  shown  in  Fig.  2.  In  this  case  the  edge  Mach  nos.  are 
heading  toward  the  sharp  cone  value,  rapidly  for  tur¬ 
bulent  flow  and  gradually  for  laminar  flow.  The  increase 
in  edge  velocity  from  the  classical  solution  will  have  a 
direct  effect  on  heat  transfer  rates  as  shown  by  the  heat¬ 
ing  equations  in  Section  4.2.  Although  the  effect  of  in¬ 
creased  velocity  significantly  increases  the  heat  transfer 
coefficient,  the  increased  heating  rate  is  moderated  by 
an  accompanying  decrease  in  the  recovery  enthalpy. 


:<>L.  Lees,  “Laminar  Heat  Transfer  Over  B.ant-Nosed  Bodies  at 
Hypersonic  Flight  Speeds,”  Jet  Propul.  26,  259-269  (1956). 


ABBOHT  can  be  run  using  either  ideal  or  equilibri¬ 
um  air  properties.  Except  when  noted  otherwise,  all 
cases  presented  herein  were  run  with  the  latter. 


15°  cone,  /?„  =  0.2b  in  ,  M  •  7.5,  TW  =  540°  R 


—Turbulent 

swallowing 

— — Laminar 
swallowing 

- No 

swallowing 

— —Sharp  cone 
valve 


Figure  2  Entropy  along  a  boundary  layer  edge. 


27E.  R,  G.  Eckert,  “Engineering  Relations  for  Friction  and  Heat 
Transfer  to  Surfaces  in  High  Velocity  Flow,”  J.  Aeronaut.  Set. 
12,  585-587  (1955). 
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15°  cone,  /?„=  0.25  in.,  M  =  7.5,  TW  =  540°  R 


S/R„ 


Rn  =  3.0  in.,  o,  =  0,  M  =  6.6,  7T  =  3260°  R, 
PT  =  1152  psia,  W  =  540° R,  QS  =  65.0  Btu/ft2-s 


—  With 
swallowing 

—Without 

swallowing 

0  NASA 
data2* 


Figure  3  Mach  no.  along  a  boundary  layer  edge. 


Figure  4  Heat  transfer  on  a  12.5*  blunted  cone. 


The  laminar  heating  distributions  for  the  old  and  new 
versions  of  ABBOHT  are  compared  to  experimental 
data  from  Ref.  28  in  Fig.  4.  Included  with  the  test  con¬ 
ditions  is  Lees’s  stagnation  point  heating  rate  (none  of 
the  experimental  data  presented  included  a  measurement 
of  stagnation  point  heating).  Both  codes  accurately  pre¬ 
dict  laminar  heating  rates  for  the  given  conditions. 
Streamline  swallowing  has  very  little  effect  on  heating 
rates  in  the  nose  region  of  this  body. 

Figure  5  shows  the  laminar  heating  distribution  for 
a  less  severe  heating  condition.29  Here  again,  the  two 
theories  fit  the  data  equally  well.  Laminar  heating 
predictions  are  within  10%  of  the  experimental  data  in 
Figs.  4  and  5. 

The  locations  most  affected  by  streamline  swallow¬ 
ing  are  downstream  of  the  nose  region  in  the  turbulent 
BL  flow  regime.  This  point  is  best  illustrated  in  Fig. 
6,  which  shows  heat  transfer  to  a  slightly  blunted  5° 
cone  at  Mach  10.  The  difference  between  the  old  and 
new  versions  of  ABBOHT  is  approximately  20%  in  the 
turbulent  flow  regime.  Results  from  both  versions  are 
compared  to  Mayne’s  method  as  well  as  experimental 
data  obtained  by  Mayne. 30  The  turbulent  heating  dis¬ 
tribution  predicted  by  ABBOHT  is  not  as  high  as 

28R.  J.  Nowak,  C.  W.  Albertson,  and  L.  R,  Hunt,  “Aerothermal 
Tests  of  a  12.5’  Con  it  Mach  6.7  for  '/arious  Reynolds  Num¬ 
bers,  Angles  of  Attack,  and  Nose  Shap*..,”  NASA  TP-2345  (Jan 
1985). 

SD.  B.  Carver,  "Heat  Transfer,  Surface  Pressure  and  Flow-Field 
Surveys  on  Conic  and  Biconic  Models  with  Boundary  Layer  Trips 
at  Mach  Number  8— Phases  IV  and  VI,”  AEDC-TSR-80-V14, 
Arnold  Engineering  Development  Center,  Tullahoma,  Tenn.  (Mar 
1980). 

WA.  W.  Mayne,  Jr.,  “Calculation  of  the  Boundary-Layer  Flo> 
the  Windward  Symmetry  Plane  of  a  Spherically  Blunted  Axisym- 
metric  Body  at  Angle  of  Attack,  Including  Streamline-Swallowing 
Effects,”  AEDC-TR-73-166,  Arnold  Engineering  Development 
Center,  Tullahoma,  Tenn.  (Oct  1973). 


Rn  »  0.5  in.,  M  ”  8, 7"7"“  1280°R,PT  =  211  psia, 
71V  *  540°  R,  QS  -  11.2  Btu/ft*-* 


Figure  5  Heat  transfer  on  a  7°  blunted  cone. 


Rn  =  0.15  in.,  M  =  10.  TT  =  2000°  R,  PT  =  1265  psia, 
71V  =  530°  R,  QS  =  64.6  Btu/ft2- s 


Figure  6  Heat  transfer  on  a  5°  blunted  cone. 
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Mayne’s  theory  at  x/R„  <  250  but  converges  with  it 
at  the  sharp  cone  value  farther  downstream.  The  tur¬ 
bulent  experimental  data  have  considerable  scatter  and 
range  from  19%  higher  to  10%  lower  than  the  new  AB- 
BOHT  predictions.  This  is  a  marked  improvement  from 
the  solution  without  swallowing,  which  ranges  from  13 
to  48%  lower  than  the  experimental  data  in  the  turbu¬ 
lent  flow  regime.  Laminar  flow  calculations  show  a 
slightly  better  correlation  to  experimental  data  for  AB- 
BOHT  results  rather  than  for  Mayne’s  theory. 

Figures  7  and  8  show  additional  comparisons  to  ■  - 
perimental  turbulent  data.  In  general,  the  curves  with 
swallowing  overpredict  heat  transfer  in  the  fully  tur¬ 
bulent  flow  regime  by  5  to  15%,  while  those  without 
swallowing  underpredict  heating  from  10  to  25%.  In 
Fig.  8,  a  BL  trip  was  used  at  an  S/Rn  of  4.  It  is  be¬ 
lieved  that  low  Reynolds  no.  turbulent  flow  (low  BL 
vorticity)  induced  significant  vorticity  interaction  in  the 


R„  =  3.0  in.,  M  =  6.6,  TT  =  3185° R,  PT  =  1230  psia, 
a  =  10,  7W  =  540°R,  QS  *  64.4Btu/ftJ-s 


- With 

swallowing 

- Without 

swallowing 

O  NASA 
data1* 


Figure  7  Heat  transfer  on  a  12.5°  blunted  cone. 


Rn=  0  5  in.,  M  =  8,  TT  =  1350°R, PT  =  857  psia, 
TW  =  540° R,  OS  =  26  5  Btu/ftz-s 


Figure  8  Heat  transfer  on  a  7°  blunted  cone. 


nose  region  of  this  body,  accounting  for  the  high  esti¬ 
mate  of  heat  transfer  at  S/Rn  <  30.  In  addition  to  im¬ 
proving  heating  rate  predicitons  by  5  to  10%,  the 
method  with  swallowing  also  produces  heating  rate  dis¬ 
tributions  more  representative  of  the  trends  in  the  data. 

To  assess  ABBOHT’s  performance  with  streamline 
swallowing  at  a  higher  Mach  no.,  two  cases  were  run 
with  ideal  air  properties  for  comparison  with  the  results 
of  the  DeJarnette  and  Zoby  methods,  obtained  from 
Ref.  31.  Additionally,  a  VSL  solution,  which  accurately 
models  variable  entropy  as  well  as  vorticity  interaction 
effects  (also  obtained  from  Ref.  31),  is  given.  Figures 
9  and  10  show  turbulent  heating  results  for  a  5°  cone 
at  Mach  15  and  angles  of  attack  of  0  and  20°,  respec- 


Rn  ■=  1.5  in.,  M »  15,  a  =  0?  TW  =  2260°R, 
alt.  “  150,000  ft,  perfect  gas 


Figure  9  Turbulent  heating  on  a  5°  cone. 


Rn  =  1.5  in.,  M=  15,  O  =  20°,  TW  =  2260°  R, 
alt  *  150,000  ft,  perfect  gas 


Figure  10  Turbulent  heating  on  a  5°  cone. 


3IE.  V.  Thompson,  E.  V.  Zoby,  K.  E.  Wurster,  and  P.  A.  Gnoffo, 
“An  Aerothermodynamic  Study  of  Slender  Conical  Vehicles,” 
AIAA  Paper  No.  87-1475,  presented  at  A1AA  22nd  Thermophys- 
ics  Conf.  (Jun  8-10,  1987). 
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tively.  In  both  cases,  ABBOHT  is  generally  higher  than 
the  DeJarnette  and  Zoby  solutions  and  below  the  VSL 
solution.  The  total  spread  in  all  four  solutions  is  ap¬ 
proximately  18%.  Usually,  a  VSL  solution  is  consid¬ 
ered  the  most  accurate  since  it  most  closely  models  the 
physics  of  fluid  heat  and  momentum  transfer.  However, 
without  experimental  data,  it  is  impossible  to  support 
this  presumption.  Nevertheless,  the  close  proximity  of 
the  ABBOHT  solution  to  the  VSL  solution  is  a  positive 
result.  It  also  is  encouraging  that  ABBOHT  accurately 
modeled  turbulent  heat  transfer  on  a  blunted  cone  with 
a  moderately  high  flow  incidence  angle  (25°)  at  high 
Mach  no.  and  altitude  conditions,  where  vorticity  in¬ 
teraction  is  moderate. 

To  appreciate  the  limits  of  the  new  ABBOHT  ver¬ 
sion,  a  case  was  run  that  included  high  vorticity  inter¬ 
action.  The  body  was  a  40°  sphere  cone  with  a  1-ft  nose 
radius.  The  flow  conditions  were  obtained  from  Ref. 
32.  For  turbulent  flow,  transition  was  assumed  to  be 
instantaneous  at  the  nose/body  tangency  point.  The 
resulting  laminar  and  turbulent  heating  distributions  are 
compared  to  Mayne’s  method  and  a  VSL  solution  from 
Ref.  32  in  Fig,  11. 

Both  ABBOHT  and  Mayne’s  turbulent  streamline 
swallowing  technique  (neither  of  which  account  for  vor¬ 
ticity  interaction)  overpredict  the  turbulent  VSL  solu¬ 
tion  by  as  much  as  60  and  40%,  respectively,  A  check 
of  the  pressure  distribution  and  entropy  swallowing  dis¬ 
tance  in  ABBOHT  revealed  them  to  be  approximately 
the  same  as  the  VSL  solution.  Sources  of  error  in  the 
ABBOHT  solution  include  the  absence  of  a  well-defined 
BL  edge  caused  by  vorticity  interaction  and  loss  of  ac- 


Rn=  1.0  ft,  M=  10,  TWA'T  =  0.1, TT  =  10161°R„ 
PT  =  571  atm  ,  QS=  130  Btu/ft2-s,  perfect  gas 


- Turbulent 

heating 

- Laminar 

heating 

—  -Mayne 


- VSL 


Figure  11  Heat  transfer  on  a  40°  blunted  cone. 


curacy  in  the  turbulent  heating  relation  at  high  incidence 
angles.  Because  of  the  limited  range  of  the  approxima¬ 
tions  employed,  ABBOHT  should  be  used  only  to  pre¬ 
dict  turbulent  heating  to  axisymmetric  bodies  with 
incidence  angles  of  less  than  25% 

Laminar  heating  predictions  were  within  3%  of  the 
VSL  solution  for  Mayne’s  method  (not  shown)  and 
within  10%  for  ABBOHT,  The  favorable  comparison 
for  laminar  flow  suggests  that  ABBOHT  can  be  used 
at  high  incidence  angles  to  predict  laminar  heating  un¬ 
der  these  conditions.  The  equivalent  cone  method  (sum¬ 
ming  the  cone  angle  and  angle  of  attack)  does  not  work 
well  for  laminar  flow31  and  should  therefore  be 
avoided. 


6.0  SUMMARY  AND  CONCLUSIONS 


A  simplified  technique  has  been  developed  to  add 
streamline  swallowing  to  an  inviscid  axisymmetric  heat 
transfer  code  named  ABBOHT,  Streamline  swallowing 
effects  are  significant  downstream  of  the  nose  region 
on  axisymmetric  blunt  bodies  where  BL  edge  entropy 

JJE.  C.  Anderson  and  D  C.  Wilcox,  “Vorticity  Interaction  Effects 
on  Blunt  Bodies,”  NASA  CR-2778  (Jan  1977). 


is  significantly  lower  than  in  the  nose  region.  To  ac¬ 
count  for  variable  BL  entropy,  a  mass  balancing  tech¬ 
nique  was  used  to  couple  the  inviscid  pressure  distribu¬ 
tion  with  mass  flow  rates  from  newly  added  BL  devel- 

,5E.  V.  Zoby  and  A.  L.  Simmonds,  “Engineering  Flowfield  Meth¬ 
od  with  Angle-of-Attack  Applications,”  J.  Spacecr,  Rockets  22, 
398  (1985). 
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opment  models.  In  addition  to  the  BL  growth  models, 
only  an  empirical  bow  shock  definition  was  required 
to  complete  the  geometry  for  the  new  procedure.  The 
new  ABBOHT  code  is  simpler  than  previous  mass  bal¬ 
ancing  techniques  because  it  does  not  solve  the  com¬ 
plete  BL  equations,  it  does  not  include  the  use  of  stream 
functions,  and  it  does  not  require  a  definition  of  the 
entire  inviscid  flow  field. 

Compressible  laminar  BL  growth  was  calculated 
directly  with  Ness’s  equations'8  for  the  BL  and  dis¬ 
placement  thicknesses  (see  subsection  3.1.1).  Compress¬ 
ible  turbulent  BL  development  was  accomplished  using 
Green’s  method,22  which  marches  the  first-order  BL 
momentum  integral  equation  together  with  the  mass  en¬ 
trainment  equation..  Both  procedures  require  iteration 
for  BL  edge  conditions. 

Comparisons  to  the  old  version  of  ABBOHT  indicate 
a  significant  increase  in  turbulent  heat  transfer  (from 
10  to  30%)  aft  of  the  nose  region.  New  ABBOHT  heat¬ 
ing  results  matched  experimental  data  almost  as  well 
as  Mayne’s  method  for  a  5°  blunted  cone  in  Mach  10 


flow.  In  general,  experimental  data  ranged  from  19% 
higher  to  10%  lower  than  ABBOHT  predictions;  vir¬ 
tually  all  of  the  scatter  was  associated  with  turbulent 
flow  cases.  Turbulent  heating  comparisons  to  DeJar- 
nette,  Zoby,  and  VSL  solutions  at  Mach  15  were  very 
good  up  to  a  flow  incidence  angle  of  25°.  However, 
at  a  larger  incidence  angle  (with  high  vorticity  interac¬ 
tion),  ABBOHT  significantly  overpredicted  heat  transfer 
rates.  Laminar  heating  was  only  slightly  higher  with 
streamline  swallowing  and  within  10%  of  the  experimen¬ 
tal  results. 

Until  further  validation  runs  and  improvements  to 
the  code  are  made,  the  Mach  15  case  with  an  incidence 
angle  of  25°  should  be  considered  the  upper  bound  of 
the  applicability  range  for  the  new  ABBOHT  code. 
Above  25°,  turbulent  heat  transfer  calculations  will  be 
high,  A  listing  of  the  ABBOHT  source  code  with  all 
subroutines  (except  the  air  property  curve  fitting  al¬ 
gorithm)  is  contained  in  Appendix  A  and  the  variable 
definitions  are  presented  in  Appendix  B. 
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APPENDIX  A 


PROGRAM  TO  COMPUTE  HEAT  TRANSFER  ON  BLUNT  BODIES 
VERSION  2.0,  JULY  1987 

ENTHALPY  BASE  OF  OUTPUT  IS  0°R,  ENTHALPY  BASE  OF 
AIRTBL  USED  IN  CODE  IS  0°R 


C  MAIN  PROGRAM 

COMMON  NOUT , NP , XTAB ( 1 00 ) , YTAB ( 1 00 ) , DSO , TO , PO , PTO , TTO , ALF , TIME , 

1  HRT(50) ,RUC(50) ,STA(50) ,XLEN(50) ,NPTS ,ZO, VO, STRAN,RN,DELC, 

2  RR.RG.LT, ISTEP,IX,SDIST(10) ,FREEM(20,10) ,EDGEMR(20,I0) , 

3  EDGEPR(20, 10) .NFIL.NUM.JSTEP , ITY, IAIR 
COMMON/ALPHAN/TITLE , HEAT , CONV , BUG, TYPE, QS, TRAN , TURB , GEO , LAM 
COMMON/LOGICAL/NURL 

CHARACTER*1  QS.NURL 
CHARACTER* 3  CONV , BUG , HEAT , LAM 
CHARACTER*4  TYPE , TRAN , TURB , GEO 
CHARACTER*7  NFORMAT(IO) 

CHARACTER*60  TITLE 

CHARACTER*80  FILES , FILE6 , FILE7 , FILE8 , FILE9 , FILE10 , FILE11 , FILE12 , 
1FILE13 , FILE14 , FILE15 , FILE16 , FILE17 , FILE18 
REAL*4  MACH.MO.LT 

DATA  NFORMAT/' FIRST' , ' SECOND* . 'THIRD' , ’FOURTH’ , 'FIFTH' , 'SIXTH' , 

1 ' SEVENTH ' , ' EIGHTH  *  , ' NINTH ' , ' TENTH' / 

C  .  NAME  INPUT  AND  OUTPUT  FILES  . 

924  F0RMAT(A80) 

925  FORMAT (3 OH  INPUT  DATA  FILE  -?  ) 

WRITE(0, 925) 

READ ( 0, 924 )FILE5 
OPEN ( 5 , FILE-FILE5 ) 

926  FORMAT (30H  OUTPUT  DATA  FILE  -?  ) 

WRITE(0 , 926) 

READ (0, 924) FILE6 
OPEN ( 6 , FILE-FILE6 ) 

928  FORMAT (30H  PLOT  OUTPUT  FILE  -?  ) 

VJRITE(0 , 928) 

READ(0, 924)FILE8 
OPEN ( 8 , FILL-FILE8 ) 

IX—  1 

C  .  READ  INPUT . . 

READ(5,910,ERR-998)  TITLE 

READ (5,914)  I AIR , ITY , IQS , ITR , IGE , NPTS , NOUT , NP 
READ(5,915)  DELC . RN , RR , RG , ALF , LT , DSO , STRAN 
IF(ITY.EQ.1)TYPE-'FLT  ' 

IF(ITY. EQ. 2 .OR. ITY. EQ  3)TYPE-'TEST’ 

IF(IQS.EQ.1)QS-'L' 

IF(IQS.EQ.2)QS-’D' 

IF(  ITR .  EQ .  1 )  TRAN- '  REYS  '■ 

I F  ( I TR .  EQ .  2 )  TRAN- '  REYT 
IF(ITR. EQ. 3)TRAN-*  PANT' 

IF(  IGE.  EQ.l)  GEO- 'CONE'' 

IF(IGE.EQ,3 )  GEO-  ’  J  LAT  '• 

IF(IGE.EQ.4;GEO-'DISK'' 

IF(IGE  EQ  5)GE0-'SPHE’- 
IF(IGE  EQ .  6 ) GEO- '  CYLI  '• 
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C . READ  NUMBER  0?  URLIM  FLOW  TABLES  TO  CREATE 

WRITE(0,1003) 

1003  FORMAT{ IX, 'CREATE  URLIM  LOCAL  FLOW  TABLES?  (Y/N)’  ) 

READ(0,1004)NURL 

1004  FORMAT (Al) 

IF(NURL.EQ, 'Y' .OR.NURL.EQ. 'y' )THEN 
WR1TE(0,929) 


C . READ  FLOW  DISTANCE  TO  EACH  STATION  (FT)  . 

929  FORMAT ( IX, 45H  NUMBER  OF  FLOW  STATIONS  (UP  TO  10)  -  ?  ) 

READ(0,930)NUM 

930  FORMAT(I5) 

WRITE(0,960) 

960  FORMAT (5X, 'LIST  SURFACE  DISTANCE  TO  EACH  STATION  (FT)  USING  FREE 
1 FORMAT  ') 

READ(0,*) (SDIST(I) ,1-l.NUM) 

NEWNUM-8  +  HUM 

DO  12  1-9.NEWNUM 

WRITE ( 0 , 939 ) NFORMAT ( I - 8 ) 

939  F0RMAT(1X,'  FILE  NAME  FOR  ',A7,'  STATION  -  ?’) 

IF(I.EQ.9)THEN 
READ(0,924)FILE9 
OPEN ( I , FILE-FILE9 ) 

ELSEIF(I . EQ, 10)THEN 

READ(0,924)FILE10 

OPEN(I.FILE-FILEIO) 

ELSEIF(I . EQ. 11)THEN 

READ(0,924)FILE11 

OPEN(I.FILE-FILEll) 

ELSEIF(I.EQ.12)THEN 
READ (0,926) FI LEI 2 
OPEN ( I , FILE-FILE12 ) 

ELSEIF(I . EQ. 13)THEN 

READ(0,924)FILE13 

OPEN(I,FILE-FILE13) 

ELSEIF(I.EQ.14)THEN 
READ(0,924)FILE14 
OPEN ( I , FILE-FILE14 ) 

ELSEIF(I.EQ.15)THEN 
READ(0,924)FILE15 
OPEN(I .FILE-FILE15) 

ELSEIF(I . EQ. 16)THEN 
READ(0,924)FILE16 
OPEN(I .FILE-FILE16) 

ELSEIF( I . EQ . 17)THEN 
READ(0,924)FILE17 
OPEN(I .FILE-FILE17) 

ELSE1F(I.EQ.18)THEN 

READ(0,924)FILE18 

OPEN(I,FILE-FILE18) 

END  IF 

12  CONTINUE 
ENDIF 


C  . READ  WALL  TEMPS  AT  NPTS  STATIONS 

IF  (NPTS . EQ  0)  GO  TO  1000 
READ(5 , 915)  (XLEN(I) , 1-1 ,NPTS) 

READ(5 , 915)  (STA(I) ,1-1, NPTS) 
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C  .  READ  INPUT  FLOW  PARAMETERS  . 

IF(ITY.EQ.1.0R.ITY.EQ.2)  GO  TO  549 
READ ( 5 , 915 , END-999 )TIME , ZO , VO , PTO , TTO 
GO  TO  551 

549  READ (5,915, END-9 99 ) TIME , ZO , VO , PO , TO 
551  JSTEP-JSTEP+1 
910  FORMAT (A60) 

914  FORMAT (101 5) 

915  FORMAT(7F10.0) 

C  .  CALL  AEROTHERMAL  HEATING  (AND  PRESSURE  DIST)  SUBROUTINE 


CALL  QDIST 
GO  TO  1000 


C  .  CONSTRUCT  URLIM  TABLES  IN  SEPARATE  FILES 

999  IF(NURL.EQ. 'Y* .OR.NURL.EQ. 'y' )THEN 
NFIL-9 

DO  11  J-l.NUM 
URITE(NFIL,901)SDIST(J) 

WRITE (NFIL, 902) 


C  .  WRITE  VALUES  IN  URLIM  TABLES  FOR  M-0,1  . 

WRITE(NFIL, 903)0., l.,l. 

WRITE(NFIL, 903)1. ,1. ,1. 

901  FORMAT (IX, 'LOCAL  FLOW  TABLE  FOR  STATION  AT  S-' , F10. 2 ,1X, 'FT' ) 

902  FORMAT (1 OX, 'MO' ,13X, 'M/MO' ,11X,’P/P0’ ) 

C . WRITE  URLIM  FLOW  TABLES  FOR  M  >  1  . 

WRITE(NFIL,903)(FREEM(I,J),EDGEMR(I,J),EDGEPR(I,J),I-1,JSTEP) 

903  FORMAT (IX, F12 .4 , 2F14.4) 

NFIb-NFIL+1 

11  CONTINUE 


ENDIF 

C  IF(IX.EQ.l)  GO  TO  900 

IF(IX.EQ.l)  GO  TO  998 
IX-IX-1 
998  END 

C . .  END  MAIN  PROGRAM  . 

C  .  AEROTHERMAL  HEATING  (AND  PRESS.  DIST.)  SUBROUTINE  - 

SUBROUTINE  QDIST 

C  /* . 

C  THIS  SUBROUTINE  (QDIST)  USES, 

C  ATM  THE  1962  STANDARD  ATMOSPHERE 

C  AIRTBL  MAIN  DRIVER  FOR  AIR  PROPrRTIES 

C  WRKHAT  THE  CURVE-FIT  AIR  PROPERTIES 

C  PIF1  LINEAR  INTERPOLATION 

C . */ 

COMMON  NOUT ,  NP , XTAB ( 100 ) ,  YTAB ( 100)  , DSO ,  TO ,  PO ,  PTO ,  TTO ,  ALF ,  TIME ,: 


1  HRT(50) ,RUC(50) ,STA(50) ,XLEN(50) , NPTS , ZO , VO , STRAN , RN , DELC , 

2  RR.RG ,LT, ISTEP, IX, SDIST(IO) , FREEM(20 ,10) ,EDGEMR(20, 10) , 

3  EDGEPR(20 , 10) , NFIL, NUM.J STEP, ITY, I AIR 
COMMON/ALPHAN/TITLE , HEAT , CONV , BUG ,  TYPE , QS , TRAN , TURB , GEO , LAM 
COMMON/LOG I CAL/NURL 

CHARACTER*!  QS.NURL 
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REAL*4  NUW,NUE,MDOT,ME2 ,M02 , MNORM , JUNK , JUNK1 , JUNK2 , JUNK4 
DIMENSION  STAB(IOO) 

C .  INITIALIZE  VARIABLES  AND  CONSTANTS . 


HREF-8465.0 
PREF-2116.224 
TGJ-50062 . 744 
TDUM-1 . 

GAM-1.40 

CPO-O.24 

RAD-57.2957795 

DTR-1 . 0/RAD 

SO-l.OE-5 

RHOO-l.O 

AO-1.0 

WO-l.O 

PW-1.0 

TS  i-l.O 

HS-1.0 

SST-1.0 

RHOST-l.O 

AST-1.0 

VST-1.0 

MMARKER-0 

MARKER-0 

I FLAG-0 

DBL1-0 . 

DBL2-0 . 

DDD-0 . 

DD1-0. 

DD2-0 . 

VE1-0 . 

VE2-0. 

TTS-1.0 

TE-1 . 0 

HU-1.0 

PR-0.72 

IFLAG1-0 

RH0S-1.0 

HE-1.0 

SW-1.0 

SS-1.0 

RH0E-1.0 

RH0U-1.0 

MUREF-1 .15312E-5 
RHREF-0. 0803707 
DIM-1.0 
1 1-1 
12-2 

14- 4 

15- 5 

PI-3.14159 

KK-0 

CONV  -  'YES' 

LAM- 'YES' 

SMAX  -  1 . 1*XLEN(NPTS) 

NN-0 

IX  -  1 

DIM  -  1.0 

RET-0.0 

DS  -  DSO 
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FSXL-0.0 

SX-0.0 

FSTL-0.0 

ST-0.0 

RES-0.0 

SR-0.0 

SL-0.0 

KKK-0 

IF  (TYPE. NE. 'TEST')  GO  TO  260 
C . INPUT  TOT  A?.  FREE  STREAM  TEST  CONDITIONS 


SSO-O. 

IF(ITY.EQ.2)G0  TO  261 
IF(IAIR. EQ. 1)THEN 

CALL  PAIR(PTO,TTO,HTO,SSO,DUl,DU2,I4) 

ELSE 

CALL  AIRTBL(PT0/144.  ,TTO,HTO,SSO,DU1,DU2,DU3,I4,N) 
ENDIF 

HO  -  HTO  -  VO**2/TGJ 

IF(IAIR.EQ.1)THEN 

CALL  PAIR(PO,TO,HO,SSO,RHOO,AO,I5) 

PO-FO/144 . 

ELSE 

CALL  AIRTBL( PO , TO , HO , SSO , RHOO , AO , DU3 , 15 ,  N) 

ENDIF 

AO-49. 014*SQRT(T0) 

MO-VO/AO 
PO— PG*144 . 

GO  TO  262 

C . INPUT  STATIC  FREE  STREAM  TEST  CONDITIONS . 

261  AO-49. 014*SQRT (TO) 

CO  TO  280 

C . INPUT  FLIGHT  CONDITIONS . 

260  CALL  ATM(ZO,PO,TO,DUl , AO, RHOO) 

280  MO-VO/AO 

C  .  CALCULATE  TOTAL  ENTHALPY  . 

HTO  -  VO**2/TCJ  +  CPO*TO 

C  .---CALCULATE  OTHER  FREE  STREAM  STATIC  AND  TOTAL  CONDITIONS 
SSO-O. 

IF(IAIR.EQ.1)THEN 

CALL  PAIR( PO , TO , HO , SSO , RHOO , AO , 14 ) 

CALL  PAIR(PTO.TTO,HTO,SSO,DUl,DU2.I5) 

ELSE 

CALL  AIRTBL( PO/144 , 0 , TO ,  HO,  SSO, RHOO, AO, WO, 14  ,N) 

CALL  AIRTBL(PTO,TTO,HTO, SSO.DUl ,DU2 ,DU3 ,15  ,N) 

PTO-PTO* 144.0 
ENDIF 


C  .  CONDITIONS  BEHIND  SHOCK 

262  IF( IAIR . EQ . 1 )THEN 

CALL  PGNS(MO , SSO, SS2 ,R2R1) 

ELSE 
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IF(IAIR.EQ.1)THEN 

CALL  PAIR(PTS,TrS,HTO,SS2,RHOS,AS,IS) 

ELSE 

CALL  AIRTBL( PTS ,TTS , HTO , SS2 . RHOS ,AS ,WS , 15 ,N) 

PTS  -  PTS*144 . 0 

ENDIF 

M02  -  MO*MO 


C  .  STAC,  FT,  VELOCITY  GRADIENT  FOR  SPHERE  --- 

BETA  -  (VO/2. /RN)*(1.+1. 56/(1.  +H02)**. 79) 

C  .  COHEN  CORRELATIONS  AT  STAG,  PT. . 

RHOMUE  -  (0 . 225*(PTS/PREF)**0 . 992)/ 

&  (1. 0-1. 0213*(1.0-(HTO/HREF)**0. 3329)) 


TW-TTS 

ssw-o. 

1F(IAIR.EQ.1)THEN 

CALL  PAIR ( PTS ,TW  HV.SStf PRH0V,AW,I4) 

ELSE 

CALL  AIRTBL(PTS/144 . ,TW,HW, SSW.RHOW, AW.WW, I4,N) 

ENDIF 

RHOMUV  -  (0. 22 5*(PTS/PREF) **0.992)/ 

A  (1. 0-1. 0213*(1.0-(HW/HREF)**0. 3329)) 

RHOMUE-RHOMUE*MUREF*RHREF 
RHOMUW-RHOHUW*MUREF*RHP.EF 
NUE-RH0MUE/RH0S**2 . 

NUW-RHOMUW/RHOW**2 , 

C  .  STAG,  PT.  BOUNDARY  LAYER  THICKNESS  . 

DDELTAS  -  6.*  0.763*  (RHOMUE/RHOMUV)**0.4 
A  *(NUW/BETA)**0.5 

C  .  VON  KARMAN  GEOMETRY  (WITH  35  DEG  NOSE  TIP)--- 

IF(GE0.NE. 'VK')  GO  TO  14 
DEL-DELC+ALF 

XT  -  RN*(1.0-SIN(DEL*DTR)) 

XTC-  RN*(1 .0-SIN(DELC*DTR) ) 

YT  -  RN* ( COS ( DEL*DTR ) ) 

YTC-  RN* ( COS ( DELC*DTR ) ) 

STN  -  RN* (90.0 -DEL) /RAD 
SRTP-YT 

CALL  SCONE (GAM , PC ,M0 , DEL , PO) 

C  . BLOG  GEOMETRY . */ 

14  IF  (GEO. NE.' BLOG')  GO  TO  300 
DELTA-ACOS ( ( RG - RR ) / ( RG - RN ) ) 

DELTA-DELTA*RAD 
DELa-DELTA  +  ALF 
RRR-(RN-RC)*COS(DEL*DTR)  +  RG 
YT-RG* ( RG - RRR ) / ( RG - RN ) +RRR - RG 
YTC-RG*(RG-RR)/(RG-RN)+RR-RG 
XT-RN - SQRT (RN**2 . -YT**2 . ) 

STN  -  2.0*RN*ASIN(SQRT(.Xl**2.+YT**2.)/2./RN) 

TH1  -  90.0-(180.0*STN/RN/PI) 

300  IF(GEO.NE.'CONE')  GO  TO  320 
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C . BLUNTED  CONE  GEOMETRY . 

XT  -RN*(1.0  -  SIN(DEL*DTR)) 

YT  -  RN  *  COS(DEL*DTR) 

YTC  -  RN  *  COS(DELC*DTR) 

XTC-RN*(1 .0-SIN(DELC*DTR) ) 

STN-RN* ( 90 . 0- DEL) /RAD 
SRTP-YT 

C . STAGNATION  POINT  VELOCITY  GRADIENTS  FROM  PERINI 

320  SS-SS2 

IF(GEO.EQ,'SPHE' .OR.GEO.EQ. 'CONE' .OR.GEO.EQ. 'ENNK' .OR. 
1  GEO.EQ. 'BLOG' .OR.GEO.EQ. ’VK' ) 

2DUEDS-(VO/2./RN)*(l. 1+1. 56/(1. +M02)**. 79) 

C  - .  CHECK  FOR  CONVERGENCE  IN  AIRTBL  . 

IF(N.LE.l)  GO  TO  330 
VRITE(6.920) 

920  F0RMAT(22H  NO  CONV,  GO  TO  NEXT  ) 

CONV  -  'NAY' 

RETURN 

C  .  BEGIN  MARCHING  . 


330  S-  SO-DS 
100  SL  -  S 

S  -  S  +  DS 

IF(GEO.EQ. 'VK' .AND.S.GE. (STN+0. 04134) )GO  TO  1505 

1504  SEFF-SEFF  +  DS 
XF-0.0 

GO  TO  1502 

C .  GEOMETRY  FOR  VON  KARMAN . 

1505  PHI1-AC0S(1.-(2.*X/LT)) 

Y1-RR*SQRT( (PHI1- (SIN(2.*PHI1) )/2 . )/PI) 

PHI2-ACOS ( 1 . - (2 . *(X+DX)/LT) ) 

Y2-RR*SQRT ( ( PHI 2 • ( S IN ( 2 . *PH 1 2 ) ) /2 . ) /PI ) 

DY-Y2-Y1 

C .  COORDINATE  TRANSFORMATION  FOR  VK  ANGLE  OF  ATTACK 

XXX1-X-Y1*TAN(ALF*DTR) 

XXX2-(X+DX) - Y2  *TAN ( ALF*DTR ) 

YYY1-(Y1/C0S(ALF+DTR) )  +  XXX1*TAN ( ALF*DTR) 
YYY2-(Y2/COS(ALF*DTR) )  +  XXX2*TAN ( ALF*DTR) 

DYYY-YYY2 -YYY1 
DXXX-XXX2 -XXXI 


C .  INCIDENCE  ANGLE  FOR  VK  AT  ANGLE  OF  ATTACK 


THb-TH 

TH-ATAN ( DYYY/DXXX )+RAD 
S-S-DS 

DS-SQRT ( DX**2+DY**2 ) 
RADIUS-Y1+DY 


S-S+DS 


X-X+DX 


Y-Y1+DY 
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1502  1F(TST.GT.1080. . AND.  IAIR.NE.  1)  GAM-1.4493-4. 535E-5*TST 
IF(TST.GT. 3600.. AND. IARI.NE.l)  GAM-1.29 


C  .  CALCULATE  LOCAL  REYNOLDS  NUMBERS  . 

1F(KKK.EQ.0)G0  TO  416 

RHOMUE  -  (0.225  *  (PE/PREF)**0.992)/ 

&  (l.C  -  1.0213  *  (1.0- (HE/HREF)**0 . 3329) ) 


RHOMUE-RHOMUE*MUREF*RHREF 
NU£-RH0MUE/RH0E**2 . 

IF(KXK.EQ.O)NUE-NUU 
416  RES-VE*S/NUE 
RET-VE*DM/NUE 

IF(TRAN.EQ. ‘REYS1 .AND.RES.LT.STEAN)GO  TO  15 
1F(TRAN.EQ,,REYT,.AND.RET.LT.STRAN)G0  TO  15 
LAM-'NO' 

1F(MARKER . LE. 2 . OR. S .  LE . STN . OR . KKK . EQ . 0)LAM-' YES ' 

GO  TO  19 

15  LAM- 'YES' 

C  . EMPERICAL  PRESSURE  DISTRIBUTIONS . */ 

C  .  SPHERICAL  PORTION  OF  VK  NOSE  . 

19  IF(GEO.NE.'VK')GO  TO  16 
IF(S.GE.STN)GO  TO  341 
TH-(S/RN)*RAD 
X  -  RN*(1 . -COS(TH*DTR)) 

PE-PTS - ( PTS - PO ; * ( SIN (TH*DTR ) **2 ) 

Y  -  RN*SIN(TH*DTR) 

SR-Y 

TH-90 . -TH 
GO  TO  101 

C .  CONE  PORTION  OF  VK  NOSE . 

341  IF(S.GE.  (STN-r. 04134) )GO  TO  342 
X  -  XT  + ( S - STN ) *COS ( DEL*DTR ) 

Y  -  YT  +(S-STN)*SIN(DEL*DTR) 

SR-Y 

RADIUS  -  YTC  +  (S-STN)*SIN(DELC*DrR) 

TH  -  DEL 

CALL  BCONE(PE,PC,DEL,RN1X|S,MO) 

GO  TO  101 

C . VK  PORTION  OF  VK  CONTOUR . 

342  IF(TH.GF. 0. 3)  GO  TO  20 
WRITE(6 , 1501) 

1501  FORMAT (  ’  BODY  SLOPE  <  0.3  DEG  -  STOP  CALCULATIONS ' ) 

RETURN 

16  IF(GEO.NE.’CONE')  GO  TO  400 

C  ---BLUNTED  CONE  PRESSURES , (BLICK  &  FRANCIS, AIAA  JR.,  V.4,#3)  */ 

IF(S.GT.STN)  GO  TO  340 
TH  -(S/RN)*RAD 
X  -  RN*(1 ,0-COS(TH*DTR) ) 

PE  -  PTS  - (PTS-PO)*(SIN(TH*DTR)**2) 

Y  -  RN*SIN(TH*DTR) 

SR-Y 
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RADIUS-YTC+(X-XTC)*TAN (DELC*DTR) 

SR-Y 

TH-DEL 

CALL  BCONE(PE , PC, DEL, RN , X , S ,M0) 
CO  TO  101 


C  . END  OF  CONE . 

C . BLOC  SPHERICAL  NOSE 


400  IF(KX.EQ.l)  GO  TO  1013 
IF(S.GT.STN)  GO  TO  410 
TH-90 . 0- ( ISO . 0*S/RN/PI ) 

X-RN*(1.0-SIN(TH*DTR) ) 

'■-RN*C0S  ( TH*DTR ) 
jR-Y 

PE-PTS- (PTS-P0)*SIN( (90 ,0-TH)*DTR)**2 
CALL  XMACH(GAM.ME,PE/PTS) 

THL  -  TH 
CO  TO  101 

C  -  NEED  CONDITIONS  EXACTLY  AT  TANGENCY  POINT,  SO  IF  WE  ALREADY 

C  .  STEPPED  PAST  IT  WE  RETURN  TO  TANG  PT  EXACTLY  . 

410  KK-1 

TH-90 . 0 - ( 180 . 0*STN/RN/PI ) 

X-RN*(1.0-S1N(TH*DTR)) 

Y-RN*COS(TH*DTR) 

RADIUS-RN*COS (TH*DTR) 

SR-Y 

PE-PTS- (PTS-PO)*SIN( (90. 0-TH)*DTR)**2 
CALL  XMACH  (GAM ,  ME ,  PE/PTS ) 

THL  -  TH 
GO  TO  101 

C . OGIVE  SURFACE . 

1013  DELT-DS*180.0/PI/RG 
TH  -  THL-DELT 
IF(TH.GE.O. 3)  GO  TO  420 
WRITE(6,919) 

919  FORMAT (  '  BODY  SLOPE  <  0.3  DEG  STOP  CALCULATIONS ' ) 

RETURN 

420  X  -  XT+2 .0*RG*SIN( (THl-TH)*DTR/2 . 0)*COS( (THl+TH)*DTR/2 .0) 

Y-YT+2 . 0*RG*SIN((THl-TH)*DTR/2.0)*SIN((THl+TH)*DTR/2 . 0) 

SR-Y 

THC-TH - ALF 
TH1C-TH1-ALF 

RADIUS-YTC+2 ,0*RG*SIN( (THlC-THC)*DTR/2 . 0)*SIN( (TH1  C+THC)*DTR/2 . 0) 
20  IF(TH.LE. 3.0)  GO  TO  430 

C  .  SHARP  CONE  PRESSURE  . 

CALL  SCONE (GAM , PC , MO , TH , PO ) 

XF-RN*(1.0-SIN(TH*DTR))+(Y-RN*C0S(TH*DTR))/TAN(TH*DTR) 

IF(GEO.EQ.'VX') 

1XF-RN*( 1 . 0- SIN (TH*DTR) )+(YYY2 -RN*COS (TH*DTR) ) /TAN(TH*DTR) 

C  .  BLUNTED  CONE  PRESSURE  FROM  SHARP  CONE  PRESSURE  . 

CALL  BCONE(PE, PC.TH.RN.XF, S ,MO) 

430  CALL  PRMY(&\M,PE,DELT,ME,PE) 

440  CALL  XMACH (GAM, ME, PE/PTS) 
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101  N-l 


c . WALL  CONDITIONS . 

CALL  PIF1 (TW , S , XLEN , NPTS , STA) 

PW  -  PE 
SSU-0. 

IF(IAIR. EQ. 1)THEN 

CALL  PAIR(PW,TW,HW, SW.RHOW.AW.I4) 

ELSE 

CALL  AIRTBL(PW/144.0,TW,HW,SSW,RH0W,AW,WU,I4  ,N) 
END  IF 


C . CHECK  FOR  TURBULENCE  --- 

IF(LAM.EQ. 'NO' )GO  TO  63 
C  .  B.L.  EDGE  CONDITIONS 


60  IF(KKK.EQ.O)  PE-PTS 
IF(IAIR. EQ. 1)THEN 
CALL  PA1R(PE,TE, HE, SS, RH0E.AE.il) 

ELSE 

CALL  AIRTBLI PE/144 . 0 . TE . HE , S  S , RHOE , AE , WE , II . N) 

ENDIF 

IF(KKK.EQ.O)  HE-HTO 
VE-MAX<1.0,TGJ*(HT0*HE)) 

VE-SQRT(VE) 

AE-49 . 014*SQRT(TE) 

ME-VE/AE 

IF(IFLAG.GE.l)  GO  TO  56 
ME2  -  ME  *  ME 
IF(S.LT.STN)  GO  TO  99 
IF(MARKER.GE.l)  GO  TO  51 

C  --  CALCULATE  EFFECTIVE  SHARP  CONE  FLOW  DISTANCE  AT  TANGENCY  PT 


C  .  SHARP  CONE  DISTANCE  TO  TANGENCY  PT  . 

SI  -  YT/S1N(TH*DTR) 

C  .  COHEN  CORRELATION  . 

RHOMUW- ( 0 . 2 2 5* ( PE/PREF ) **0 . 992 ) / 

&  (1. 0-1. 0213*(1.0-(HW/HREF)**0. 3329)) 

RHOKUW-RHOMUW*MUREF*RHREF 
NUW-RHOMUW/RHOW**  2 

C . PSEUDO  REYNOLDS  NO,  . 

SQRE  -  SQRT(VE*S1/NUW) 

C . BL  THICKNESS  AT  TANG,  PT  -  SHARP  CONE.  LAMINAR  -  N.  NCSS- 

DDS  -  (TE*S1/1 . 732/SQRE/TV)*(3 .0*TW/TE  +  2.01* 

6  (1 .+(GAM- 1 . )*ME2/2 . -TW/TE)  -  1 .62*(GAM-1. )*ME2/2 . ) 

C  .  ADJUST  FLOW  DISTANCE  FOR  BLUNTNESS  EFFECT  . 


FF-  1. 

DBL  -  DDELTAS  +  FF  *  DDS 
VENUW-SQRT(VE/NUW) 

SEFF  -  (1 . 732*VENUW*TV*DBL/TE)/(3 . 0*TW/TE  +  2.01* 


30 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

LAUREL.  MARYLAND 


c .  CALC  BL  MASS  FLOW  FOR  LAMINAR  CASE  (LBM/S) 

51  MARKER  -  MARKER  +  1 

RH0MUW-(0 . 225*(PE/PREF)**0 .992)/ 

$  (1. 0-1. 0213*(1.0-(HW/HREF)**0. 3329)) 

RHOMUW-RHOMUW*MUREF*RHREF 
NUW-RH0MUW/RH0W**2 . 

56  ME2  -  ME  *  ME 

C .  B.L.;  MASS  FLOW  RATE . 


SQRE-SQRT (VE*SEFF/NUW) 


MOOT  -  2 . *PI*RADIUS*VE*RHOE*SEFF*TE/l . 732/SQRE/ 

&  TW  *  (2.007*TW/TE  +  2 . 01*(1.+(GAM-1 . )*ME2/2. -TW/TE) 

6.  -2.003  *  (GAM-1. )*ME2/2.) 

C . RADIUS  OF  FREE  STREAM  TO  MATCH  MDOT . 

57  YS  -  SQRT(MDOT/PI/RHOO/VO) 

C  .  SHARP  CONE  ATTACHED  SHOCK  ANGLE  . 

ETA  -  SQRT(l./MO**2.+(GAM+l.)*(SIN(TH*DTR))**2./2.) 

ETA  -  ASIN(ETA)*RAD 

C  .  BOW  SHOCK  ANGLE  (GAMMA)  -  B1LL1G  . 

RM1  -  0. 54/(MO-l . )**1.2 
RC  -  1 . 143*RN*EXP (RM1 ) 


TANGAM  -  RC/YS  *  SQRT(1.+(YS*YS*(TAN(ETA*DTR)**2 . )/RC**2 . )) 
GAMMA-RAD* AT AN ( TANGAM ) 

. PASS  THRU  BOW  SHOCK . 

VNORM  -  V0*S1N(GAMMA*DTR) 

VTAN  -  VO*COS (GAMMA*DTR) 

MNORM-VNORM/AO 

IF(IAIR. EQ. 1)THEN 

CALL  PGNS(MNORM,SSO,SS2,R2Rl) 

ELSE 

CALL  SHOCK(MNORM , VNORM , VTAN . RHOO , PO/144 . , HO . V2 , SS2 , R2R1 . P2 ) 

ENDIF 

C  .  CHECK  FOR  TURBULENCE  . 

IFOAM.EQ.  ’YES’)  GO  TO  64 
IF(IA1R.EQ.1)THEN 

CALL  PAIR(PE,TE. HE. SS2, RHOE.AE.il) 

ELSE 

CALL  AIRTBL(PE/144  ,TE,HE , SS2 ,RHOE, AE.WE, II ,N) 

ENDIF 

VE-SQRT ( TG J  * ( HTO - HE ) ) 

DVE-VE-VE2 
AE-49 ,014*SQRT(TE) 

ME-VE/AE 

C  . . COMPARE  ENTROPY  DOWNSTREAM  OF  SHOCK  W/PREVIOUS  BL  ENTROPY 

64  IF(ABS(SS-SS2).GT.0.001)THEN 
IFLAG-IFLAG+1 
1F(IFLAG.LE  10)G0  TO  49 
IFLAG1-IFLAG1+1 
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WRITE(6,600)S 
GO  TO  150 

600  FORMAT (//,  'STREAMLINE  SHALLOWING  FAILED'  ,F5. 4,  'FT  FROM 
&  NOSE  DUE  TO  AN  ENTROPY  DISCONTINUITY' ) 

C  .  SET  ENTROPY  PREVIOUS  VALUE  EQUAL  TO  CURRENT  VALUE  - 

49  SS-SS2 

IF(LAM.EQ.'NO')GO  TO  74 

GO  TO  60 

ELSE 

48  SS-SS2 
VE1-VE2 
VE2-VE 
IFLAC-0 
IFLAG1-0 

IF(LAM.EQ. 'NO* )G0  TO  99 

C . LAMINAR  BL  THICKNESS  -  N.  NESS  . 

DBL-(SEFF*TE/1 . 732/SQRE/TW)*(3 . 0*TW/TE  +  2.01*(I.+(GAM-1. ) 

&  *ME2/2 . -TW/TE) - 1 . 6 2* (GAM- 1 . )*ME2/2 . ) 

DBL1-DBL2 

DBL2-DBL 

DDD-(SEFF*TE/1 . 732/SQRE/TW)*(0 . 993*TW/TE+0 . 383*(GAH-1.  )*ME2/2 . ) 

DD1-DD2 

DD2-DDD 

DM-0 . 383*SEFF*TE/1 . 732/SQRE/TV 

GO  TO  99 

ENDIF 

C  .  TURBULENT  STREAMLINE  SWALLOWING  . 

63  MARKER-MARKER* 1 

IF(NMARKER.EQ.l)  GO  TO  61 
MMARKER-1 


. CALC  INITIAL  CONDITIONS  . 

- INITIAL  MASS  FLOW  THICKNESS  PROM  LAMINAR  MOOT 

BIGD  -  MDOT/2. /PI/RADIUS /RHOE/VE 

C  .  DISPLACEMENT  THICKNESS  . 

DD  -  DBL-BIGD 

C . MOM  THICKNESS  . . 

DM  -  RET*MUE/RHOE/VE 

C . INITIAL  ENTRAINMENT  RATE . . 

F-  ((DBL-DBL1) - (DDD-DD1) )/DS 

C . INITIAL  VELOCITY  CHANGE . 

DVE  -  VE-VE1 

C . INITIAL  SHAPE  PARAMETER . . 


H  -  DD/DM 
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C . STEP  ENTRAINMENT  EQ . 

61  DBIGD  -  F*DS  +  (ME2 - 1 . ) *BIGD*DVE/VE 

C . STEP  MOMENTUM  INTEGRAL  EQ . . 

DDM  -  CF*DS/2 .  -  (H+2.-ME2)  *  DM  *DVE/VE 

C . NEW  MASS  FLOW  THICKNESS . 

BIGD  -  BIGD  +  DBIGD/2.25+*. 2 

C  . NEW  MOMENTUM  THICKNESS  . 

DM  -  DM  +  DDM/2 . 25**. 2 

C  .  MODIFIED  COMPRESSIBLE  SHAPE  PARAMETER 


HI  -  BIGD/DM 


C . TRANSFORMED  COMPRESSIBLE  SHAPE  PARAMETER 

HBAR  -  1,  +  (0. 9/(Hl*3 . 3))**0.75 
C  .  KINEMATIC  (INCOMP)  SHAPE  PARAMETER 


TERM1  -  (HBAR+1. )* (HBAR- 1 .  )**2  .• 

TERM2  -  TERM1/HBAR/(3.*HBAR-1.)/(2.*HBAR-1.) 
TERM3  -  R*ME2*TERM2/5.: 

TERM 4  -  HBAR* (HBAR+1 . ) 

TERM5  -  TERM4/(3.*HBAR-1.)/(2.*HBAR-1.) 

TERM6  -  1.-TERM5 

TERM7  -  R  *  ME2  *  TERM6/5 . 

TERM8  -  1,  +  TERM7 
TERM?  -  TERM3/TERM8 
HK  -  HBAR  *  (1.+TERM9) 


C . KINEMATIC  MODIFIED  SHAPE  PARAMETER 

H1K  -  3.4  +  1 . 87/(HK-0. 5)**3 . 8 

C . NEW  ENTRAINMENT  RATE . 

F  -  0.0306  *  (H1K-3 ,0)**(- . 653) 

C . NEW  SHAPE  PARAMETER . . 

74  HR  -  HE  +  R*VE**2/TGJ 


IF(IAIR.EQ. 1)THEN 

CALL  PAIR(PW,TR,HR, JUNK.JUNK1 , JUNK2 ,12) 

ELSE 

CALL  A1RTBL(PW/144. ,TR,HR, JUNK, JUNK1 , JUNK2 , JUNK4 , 12 ,N) 
ENDIF 

H-(TW/TE)*HBAR+TR/TE- 1 . 

C  .  NEW  DISP.  THICKNESS  . 


DD-H+DM 


C 


NEW  BL  THICKNESS 


DBWU+DM+DD 
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C  .  INCOMPRESSIBLE  TRANSFORMED  REYNOLDS  NO  --- 

RETINC  -  (MUE/MST)  *  RET 

C . NEW  SKIN  FRICTION  FROM  INCOMP  CORRELATION 

CFINC  -  0. 246*EXP(-1. 56*HBAR)*RETINC**(-0.268) 

CF  -  (TE/TST)  *  CFINC 

C . MASS  FLOW  RATE . 

MDOT  -  2 . *PI*RADIUS*VE*RHOE*BIGD 
GO  TO  57 


C . REC  FACTOR . 

99  R  -  SQRT(PR) 

IF  (LAM. NE. ’YES’)  R-PR**( 1.0/3 . 0) 

C  . REC  ENTH . 

HR  -  HE  +  R*V£**2/TGJ 

C . ECKERT'S  REFERENCE  ENTHALPY . 

HS  -  HE  +  0. 5*(HW-HE)  +  0.22*(HR-HE) 
C . REFERENCE  CONDITIONS 


SST-0. 


IF(IAIR.EQ.1)THEN 

CALL  PAIRAW.TST.HS.SST.RHOST.AST.H) 

ELSE 

CALL  AIRTBL( PW/144 . 0 , TST , HS , SST , RHOST , AST , WST , 12  ,N) 
ENDIF 


C .  COHEN  CORRELATIONS . 

RMER  -(0. 225*(PE/PREF)**0. 992)/ 

1  (1.0  -  1.0213*(1.0  ■  (<HE+  0. 0) /HREF) **0. 3329) ) 

RMSR  -(0.225*(PE/PREF)**0.992)/ 

1  (1.0  -  1.0213*(1.0  •  ((HS+  0.0)/HREF)**0.3329)) 

RMS  -  RMSR*MUREF*RHREF 
RME  -  RMER*MUREF*RHREF 
MST  -  RMS/RHOST 
MUE  -  RME/RHOE 
IFOAM.EQ.  'NO'  )G0  TO  12 
IF(KKK.NE.O)  GO  TO  520 
RMES  -  RME 
RMSS  -  RMS 


C  .  STAG.  PT.  VELOCITY  GRADIENT  . 

DUEDSS-(VO/2 . /RN)*(1 . HI . 56/( 1 .+M02)** . 79) 

RMWR  -(0 . 225*(PE/PREF)**0 . 992 )/ 

1  (1.0  -  1.0213*(1.0  -  ( (HV+  0.0)/HREF)**0.3329)} 
RMWS  -  RMWR*MUREF*RHREF 

C . LEES  STAG  PT  COLD  WALL  (540  R)  Q - 

QSTAGL-0 . 88*SQRT(RMSS*DUEDSS)*(HTO-128 .7) 
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STAGD-17600 . 0*SQRT(RHOO/0 . 07647)*(VO/26000 . 0)**3 . 15/SQRT(RN) 
1  *((HTO-128.7)/HTO) 

QSTAGD-QSTAGD* ( 2 . 0**< (DIM-1 .0)/2 . 0) )  *SQRT(DUEDS/DUEDSS) 

C .  FAY, RIDDELL  STAG  PT  COLD  WALL  Q . 

QSTAGF-0.94*(RMWS**0.1)*(RMES**0.4)*SQRT(DUEDSS)*(HTO-128.7) 
QSTAGF-QSTAGF*(2 ,0**( (DIM-1 .0)/2 . 0) )*SQRT(DUEDS/DUEDSS) 

C . STAG.  PT.,  HEAT  TRANSFER  COEFF,  . 


RUCH1  -  QSTAGL/(HTO- 128 . 7 ) 

RUCHL-RUCH1 

HRL-  HR 

MRL-0 . 0 

PEL-PE 

GO  TO  102 

520  RES  -  RHOE*VE*S/MUE 


C  FSX1-RHOE*VE*(MUE**0 . 25)*(SR**(1 . 25*DIH) ) 

C  SX1  -  SX1+ (FSX1+FSX1L) *DS/2 . 0 

C  FSX1L-FSX1 

FSX  -  RMS*VE*(SR**(2 .0*DIM) ) 

SX  -  SX+(FSX+FSXL)*DS/2 . 0 
FSXL-FSX 

C  FST-FSX*RME/RMS 

C  ST  -  ST+ ( FST+FSTL) *DS/2 . 0 

C  FSTL-FST 

C . LAMINAR  FLOW . 

C . LEES  EQUATION  FOR  HT,  TRANSFER  COEFF. 


RUCH1-QSTAGL*RMS*VE*SR**DIM/SQRT(RMSS*DUEDS  *2 . 0*SX)/ 
1  (2.0**(DIM/2.0)  )/(HTO-128.7) 

GO  TO  102 


C  . TURBULENT  FLOW . 

C .  SHARP  CONE  HT,  XFER  COEFF.  NASA  TECH  PAPER  2345  JAN.  1985  -- 

12  RUCH1  -  0.0296*2. 25**. 2*RH0ST*VE/PR**0. 667 
1  / ( RHOST*VE*S/MST ) **0 . 20 

C  .  COLD  (ISOTHERMAL)  WALL  AND  PSEUDO -HOT  WALL  HT,  XFER  . 

102  QHW-RUCH 1*(HR-HW) 

QCW-RUCH1*(HR-128.7) 

C . PRINTOUTS . 


IF(KKK.NE.O)  GO  TO  540 
WRITE(6 ,900)TITLE 
WRITE (0 , 900)TITLE 

900  FORMAT ( IX, A6 2/) 

IF(  GEO.EQ. ’BLOG')WRITE(6,501) 

901  FORMAT (3X, 'BLUNTED  OGIVE, AXISY-FLOW’  ) 

IF(  GEO.EQ. ’CONE’ )WRITE( 6, 904) 

904  FORMAT(3X, ’AXISYMMETRIC  FLOW(CONE)’  ) 

IF(  GEO.EQ. ’VK’)WRITE(6, 1506) 

1506  FORMAT (3X, 'VON  KARMAN ,AXISYM-FLOW'  ) 

WRXTE(6,909) 

909  FORMAT (3X, ’LEES  STAG  POINT  Q*  ) 

WRITE(6,912) 

912  FORMAT ( 3X, ’TURB  Q  VIA  MOM.  INTEGRAL  EQ.  USING  1/7TH  POWER  FOR  VEL. 
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C  HOX-HO-128.7 

C  .  WRITE  TO  OUTPUT  FILE  AND  SCREEN  . 

WRITE(6,1014) 

WRITE  ( 6 , 91 3  )  TIME ,  ALF ,  DELC ,  RN ,  RR ,  RG ,  XT ,  YT 
WRITE  ( 0 , 913  )TIME ,  ALF ,  DELC ,  RN ,  RR ,  RG ,  XT ,  YT 
IF(TYPE.EQ.’TEST')GO  TO  817 

WRITE(6,1009)ZO,VO,MO,SSO,TTO,TO,PTO/2116.224,PO/2116.224,HTO,HO, 

&TW 

WRITE(0,1009)ZO,VO,MO,SS0,TTO1TO,PTO/2116.224,PO/2116.224,HTO,HO, 

&TW 

GO  TO  818 

817  WRITE  (6 , 1010)MO ,  VO ,  TTO ,  TO ,  PTO/2116 . 224 ,  PO/2116 . 224 ,  HTO ,  HO ,  TW,  SSO 
WRITE(0,1010)MO,V01TTO>TO,PTO/2116.224,PO/2116.224,HTO,HO,TW,SSO 

818  WRITE(6 . 1034) PTS/2116 . 224 , SS2 , TTS , R2R1 .P2/2116 . 224 
WRITE (0 , 1034) PTS/2116 . 224 , SS2 ,TTS ,R2R1 .P2/2116 . 224 

1014  FORMAT (IX,' ENTHALPY  -  0  AT  0  DEG  R'/ 

1  IX,  'COLD  WALL  Q  BASED  ON  WALL  TEMP  OF  540  DEG  R’/ 

2  IX, 'HOT  WALL  Q  BASED  ON  SUPPLIED  WALL  TEMP’//) 

913  FORMAT ( 

3  IX, 'TIME  (SEC)  ' ,F12 .4/ 

4  lX.’RADOME  GEOMETRY:’/ 

5  IX, '  ANCLE  OF  ATTACK  -  '.F10.3/ 

6  IX,'  CONE  ANGLE  (DEG)  -  ’.F10.3/ 

7  IX,'  NOSE  RADIUS  (FT)  -  ’.F12.4/ 

8  IX,'  BASE  RADIUS  AT  TANG.-  FT.  (FT)  -  '.F12.4/ 

9  IX,'  OGIVE  RADIUS  (FT)  -  '.F12.4/ 

1  IX,'  COORDINATES  OF  NOSE/BODY  TANG.:  PT.,  (FT)  ARE  X  -'.F5.5/ 

2  IX,'  Y-',F5.5/) 

1009  FORMAT ( 

1  IX .'FREE  STREAM  FLIGHT  CONDITIONS:’/ 

2  4X, '  ALTITUDE  (FT)  '  ,F12. 4, 2X, 'VELOCITY  (FT/S)  \F12.4/ 

3  4X, '  MACH  NO.  ' ,F12. 4, 2X, 'ENTROPY  (BTU/LBM-R)  ’.F12.5/ 

4  4X,  *  TOT.  TEMP  (R)  ’  ,F12 . 4 , 2X,  ’  STATIC  TEMP  (R)  '.F12.4/ 

54X,'  TOT.  PRESS  (ATM)  ' ,F12. 4, 2X, 'STATIC  PRESS  (ATM)  '.F12.4/ 

6  4X, '  TOT.  ENTH  (BTU/LBM) ' ,F12. 4 , 2X, 'STATIC  ENTH  (BTU/LBM) ’ ,F12 .4/ 

7  4X, '  WALL  TEMP  (R)  ’.F12.4/) 

1010  FORMAT ( 

1  IX, 'FREE  STREAM  TUNNEL  CONDITIONS:’/ 

2  4X,  ’  MACH  NO.  ' ,F12. 4, 2X, 'VELOCITY  (FT/S)  ’.F12.4/ 

3  4X, '  TOT.  TEMP  (R)  ’ ,F12. 4, 2X, 'STATIC  TEMP  (R)  ’.F12.4/ 

4  4X, '  TOT.  PRESS  (ATM)  ' ,F12 .4 , 2X, 'STATIC  PRESS  (ATM)  ’.F12.4/ 

5  4X, '  TOT.  ENTH(BTU/LBM)  ' ,F12. 4, 2X, 'STATIC  ENTH  (BTU/LBM) F12 .4/ 

6  4X, '  WALL  TEMP  (R)  ' ,F12 .4, 2X, ' ENTROPY  (BTU/LBM-R)  '.F12.5/ 

7  ) 

1034  FORMAT ( 

1  IX, ’STAGNATION  POINT  CONDITIONS:'/ 

2  4X, ’  TOT.  PRESS  (ATM)  F12 . 4 , 2X, ' ENTROPY  (BTU/LBM-R)  ’.F12.5/ 

4  4X, '  TOT.  TEMP  (R)  '.F12.4/ 

4  4X, '  DENSITY  RATIO  ACROSS  NORMAL  SHOCK  (RH02/RH01 ) ' , 8X , '  '.F12.4/ 

5  4X, '  STATIC  PRESS  BEHIND  NORMAL  SHOCK  (ATM)',8X,’  '.F12.4/ 

6) 

WRITE(6 ,611) 

WRITE(0,611) 

WRITE(6,914)QSTAGL,QSTAGF,QSTAGD 

WRITE(0,914)QSTAGL,QSTAGF,QSTAGD 

WRITE(0,926) 

611  FORMAT(2X, 'COLD  WALL  QSTAG  FOR  LEES  FAY-RIDDELL  DETRA- ’ 

1 , ' KEMP-RIDDELL(BTU/FT2 - S) ' ) 

914  FORMAT(21X,F8 .2 ,4X,F8 . 2 , 10X.F8. 2/) 
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FREEM(JSTEP,IST£P)  -  MO 
EDGEMR(JSTEP, ISTEP)  -  ME/MO 
EDGEPR ( J STEP , ISTEP)  -  PE/PO 
XSTEP  -  ISTEP  +  1 
ENDIF 

C  -  PRINT  COLUMN  HEADING  IN  OUTPUT  FILE  EVERY  40TH  WRITE  STEP- 

724  IF(KKK.EQ.O.OR.MOD<KKK,NOUT*40).EQ.O)CO  TO  560 
GO  TO  570 

560  IF(KKK.NE.O  )  WRITE(6,918) 

918  FORMAT(lHl) 

WRITE(6,915) 

915  FORMAT  (  54H  S/RN  X/RN  Y/RN  RAD/RN  TH  PE/PO  ME/MO  VE  , 

1  48H  SS  HE/HTO  RES  RET  DBL  RUCH  . 

2  31H  HR  QCW  QHW  LAM  ) 

926  FORMAT (  45H  S/RN  PE/PO  ME  TE/TTS  RUCH  HR  , 

1  17H  QCW/QS  LAM  ) 

.  WRITE  TO  OUTPUT  FILE  WHEN  REMAINDER (NOUT/KKK)  -  0  . 

570  IF(MOD(KKK,NOUT) .EQ.O.OR.KKK.EQ.O)  GO  TO  580 
CO  TO  590 
C  580  HEX-HE-128.7 
C  HRX-HR-128.7 

580  CP-64. 4*(PE-PO)/(RHOO*VO**2.) 

ST-QHW/RHOE/VE/(HR-HW) 

IF(GEO.EQ.'VK' .AND.S.GE. (STN+. 04134) )THEN 
WRITE(6,916) 

1S/RN ,  XXX2/RN ,  YYY2/RN .  RADIUS/RN .  TH ,  PE/PO ,  ME/MO ,  VE ,  SS ,  HE/HTO , 

2RES . RET . DBL , RUCH1 , HR , QCW , QHW . ST , LAM 
ELSE 

WRITE<6,916) 

1S/RN, X/RN , Y/RN . RADIUS/RN . TH , PE/PO , ME/MO , VE . SS . HE/HTO , 

2RES . RET . DBL . RUCH1 . HR . QCW , QHW . ST , LAM 
ENDIF 

916  F0RMAT(1X,2F7 . 2 ,2F6 . 2 ,F6, 1, 2F7. 3 ,F7 . 0, F7 .4 ,F7 . 3,2(1PE9. 2)  ,E10. 3 , 

1  0PF8 . 5,F8 . 1 . 2F9 . 3.E12.4, 2X,A3) 

C  - TURBULENT  DEBUG  WRITE  . 

C  WRITE(6 , 955) 

C  955  FORMAT (  52H  S/RN  PE/PO  VE  BIGD  HBAR  MDOT, 

C  1  49H  DM  DBL  F  H 

C  2  12H  HI  ) 

C  WRITE(6 , 956) 

C  1S/RN, PE/PO, VE, BIGD, HBAR, MDOT, DM, DBL, F,H,  HI 
C  956  F0RMAT(1X,2F6. 2 ,F7 .0, 8E12 .3) 

C - FOR  PLOT  FILE . 

IF(GEO.EQ.'VK’ .AND . S . GE. <STN4 . 04134) )  THEN 
WRITE(8,916) 

1S/RN , XXX2/RN , YYY2/RN , RADIUS/RN , TH , PE/PO , ME/MO , VE , SS , HE/HTO , 

2RES , RET , DBL , RUCH1 , HR , QCW , QHW , ST , LAM 
ELSE 

WRITE(8 , 916) 

1S/RN , X/RN , Y/RN , RADIUS/RN , TH , PE/PO , ME/MO , VE , SS , HE/HTO , 

2RES , RET , DBL , RUCH1 . HR , QCW , QHW , ST , LAM 
ENDIF 


C .  FOR  PRINT  TO  SCREEN 
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C . INTERPOLATION  FOR  RUCH  &  HR 

550  IF<MOD(KKK,NOUT).EQ.O)GO  TO  551 
C  HEX-HE-128.7 

C  HRX-HR-128.7 


551  C  -  (XLEN(IX)-SL)/(S-SL) 

PEX-PELfC*(PE-PEL) 

RUC(IX)  -  RUCHL  +C*(RUCH1-RUCHL) 

HRT(IX)  -  HRL  +  C*(HR-HRL) 

MRX-MRL+C*(ME-MRL) 

C  .  WRITE  TO  SCREEN  . 

WRITE(0,930) 

1 S/RN , PE/PO ,  ME . TE/TTS , RUCH1 , HR , QCW/QSTAGL . LAM 
930  FORMAT(lX,F8.2,F9.2,F7.2,F6.2,F9.5,F7.1,F6.2,4X,A3) 

C . STEP  STATION  FOR  WALL  TEMP . 

IX  -  IX+1 
105  KKK  -  KKK+1 
RUCHL  -  RUCH1 
PEL-PE 
HRL  -  HR 
MRL-ME 

IF(S  +  DS.LT.SMAX)CO  TO  100 
DS  -  SMAX-S 

IF(DS.GT.O.O)  GO  TO  100 
LAM  -  'YES’ 

150  RETURN 
END 


SUBROUTINE  SCONE(GAM,PC,MI ,DEL,PO) 

REAL*4  MI2 ,MI4 ,MI6 ,MI 
DTR-1 . 0/57 . 3 
C  GAM  -  1.4 

MI2  -  MI*MI 

C  /*  SHARP  CONE  PRESSURE  (SIMON  A  WALTER  ,  AIAA  JR.,  V.l,#7  )  */ 

SD2  -  SIN(DEL*DTR)**2 
MI4  -  HI2*MI2 
Ml 6  -  MI4*MI2 
GP7  -  CAM+7.0 

FI  -  GP7/4.0  -  (GAM-1. 0)**2/16.0  +  6.0/MI6  +• 

1  (MI2-1 .0)/MI4/SIN(DEL*DTR) 

CPC-F1*SD2 

PCO  -  1.0  +  GAM+MI 2*CPC/2 . 0 

PC  -  PCO*PO 

RETURN 

END 

SUBROUTINE  BCONE( PE , PC , DEL, RN , X , S , MO) 

C  /*  BLUNTED  CONE  PRESSURES , (BLICK  A  FRANCIS, AIAA  JR.,  V.4,#3)  */ 

REAL*4  MO, NX 
PI-3.14 

SP  -  X*(DEL*PI/180.0)**2/RN 
IF  (SP.GT.0.106)G0  TO  40 
NX-  -0. 30S*(  LOG(SP)  +  2.24) 

PE  -(0 , 929* (MO/2 , 8)**NX)*PC 
GO  TO  100 
40  CONTINUE 

IF(SP.LE. 0.106)  GO  TO  90 
IF(SP.GE.1.06)  GO  TO  90 

NX— 0.116*SQRT(MO-2.8)*SIN(  (LOG(SP)+2.24)/0.731) 

PE  -  PC*(SP**0.032)*EXP(NX) 

GO  TO  100 
90  PE  -  PC 
100  RETURN 
END 
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SUBROUTINE  PRMY(GAM,PE,DTX,M,P) 

C  /*  PRANDTL  MEYER  TURNING  P-F(DTHETA)  .NACA1135  ,EQ(174) 

C  FOR  SMALL  TURNING  ANGLES  */ 

REAL*4  M,M2,M4,M6,M2M1 
C  GAM-1.4 
M2-M*M 
M4-M2+M2 
M6-M2+M4 
GP1-GAM+1 . 0 
M2M1-M2 -1.0 
DT-DTX+3 . 14/180  • 

PR-1 . 0-GAM*M2*DT, SQRT(M2M1) 

1  +GAM*M2*(GPl*M4-4 . 0*M2Ml)*DT*DT/(4 . 0*M2M1**2 ) 

2  - (GAM*M2*DT**3/2 . 0/M2Ml**3 ,5)* 

3  (GPl*M4**2/6 . 0  - ( 5 . 0+7 . 0*GAM- 2 . 0*GAM**2 )*M6/6 . 0 

4  +5 . 0*GPl*M4/3 . 0 -  2 . 0*M2+4 .0/3.0) 

PE-PR*P 

RETURN 

END 

SUBROUTINE  XMACH(GAM,M,PPT) 

C  /*  MACH  FUNCTION,  M-F(P/PTOTAL)  NACA1135,,  EQ(44)  */ 

REAL*4  M 
C  GAM-1.4 

iMl-GAM-1.0 

M-SQRT ( ( PPT** ( -GM1/GAM) -1.0)+2. 0/GM1 ) 

RETURN 

CO  100,1  -  2 ,N 
IF(X.LE.XLIST(I))  GO  TO  200 
100  CONTINUE 
I  -  N 

200  F-FLIST(I-1)+(X-XLIST(I-1))*(FLIST(I) -FLIST(I-l))/ 

1  (XLIST(I) -XLIST(I-l)) 

RETURN 

END 

SUBROUTINE  PGNS(M1,SS1,SS2,RH0R1) 

C  -  PERFECT  GAS  •  NORMAL  SHOCK  SUBROUTINE 

REAL*4  Ml. Ml 2 
M12-M1*M1 

P20P1-  1.  +  2. 8/2. 4  *  (Ml 2  -  1.) 

R20R1  -  2.4  *  M12  /(2.+.4*M12) 

T20T1-P20P1/R20R1 

CP-.24 

R-. 068554 

SS2-SS1+CP*AL0G(T20T1)-R*AL0G(F20P1) 

RETURN 

END 

SUBROUTINE  PAIR(P,T,H,SS ,RHO,A, I) 

C  .  PERFECT  GAS  AIR  RELATIONS  . 

PREF-2116 . 224 
TREF-537, 

SREF-0. 59945 
CP- . 24 
Rl-53.34 
R2-. 068554 

GOTO( 10, 20,40,40 , 50)1 

10  T-TREF+EXP ( ( SS • SREF+R2+ALOG ( P/PREF) ) /CP) 

H-CP*T 
GO  TO  60 
20  T-H/CP 

SS-CP+ALOC (T/TREF) - R2+ALOG ( P/PREF)+SREF 
GO  TO  60 

40  SS-SREF+CP*ALOG(T/TREF) -R2*ALOG( P/PREF) 

H-CP*T 
GO  TO  60 
50  T-H/CP 

P-PREF*EXP( (CP*ALOG(T/TREF) -SS+SREF)/R2) 
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60  RHO-P/R1/T 

A-4C . 014*SQRT(T) 

RETURN 

END 


C . REAL  GAS  SUBROUTINE  TO  PASS  THROUGH  A  NORMAL  OR  OBLIQUE 

C  .  SHOCK  GIVEN  PROPERTIES  ON  UPSTREAM  SIDE 

SUBROUTINE  SHOCK (MNORM , VNORM , VTAN , RHOl , PI , HI , V2 , SS2 , R2R1 , P2P) 
REAL  MNORM , VNORM , VTAN , RHOl, PI, HI 
R2R1-6 . *MN0RM**2 . / (MNORM**2 . +5 . ) 

R2-R2R1*RH01 

ERR-100. 

DO  1C  i-1,50 

P2-P1+RH01*VN0RM**2 .*(1. -l./R2Rl)*2 . 1584E-4 

H2-H1+VN0RM**2 . /2 . *(1 . - 1 . /R2R1**2 . )/25051 . 6 

PP-P2 

HP-H2 

N-0 

SS2-0. 

12-2 

CALL  AIRTBL(PP,T2,HP,SS2,RHOP,A2,U2,I2,N) 

RH02-RH0P 

ERR-ABS ((RH02-R2)/RH02) 
lF(ERR.LT.1.0E-4)  GO  TO  11 
R2-RH02 

R2R1-RH02/RH01 

10  CONTINUE 

11  CONTINUE 
V2N-VNORM/R2R1 

V2-SQRT ( V2N**2 . +VTAN**2 . ) 

P2P-P2*144. 

RETURN 

END 


SUBROUTINE  SNUGS (PC, TR, I) 

DIMENSION  PC(3) ,TR(3) , A(2) 

PC(3J-  PC(1)  +  (1.0-TR(1))*(PC(2)-PC(1))/(TR(2)-TR(1)) 
A(l)  -  ABS(TR(1)-1.0) 

A(2)  -  ABS(TR(2) -1.0) 

TMAX  -  MAX  (A(l) ,A(2)) 

DO  100  I  -  1,2 
IF  (TMAX.EQ.A(I))GO  TO  9 
100  CONTINUE 
9  PC(I)  -  PC<3) 

RETURN 

END 
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APPENDIX  B 

ABBOHT  VARIABLE  DEFINITIONS 


AE 

ALF 

ALT 

AO 

AS 

AST 

AW 

BETA 

BIGD 

CF 

CFINC 

CONV 

CP 

CPO 

DBIGD 

DBL 

DBL1.2 

DD 

DD1,2 

DDD 

DDELTAS 

DDM 

DDS 

r»ci 

L/i-ii-i 

DELC 

DELT 

DELTA 

DIM 

DM 

DS 

DTR 

DU1 

DU2 

DU3 

DUEDS 

DUEDSS 

DVE 

DXXX 

DYYY 

EDGEMR 

EDGEPR 

ETA 

F 

FF 

FREEM 


speed  of  sound  at  BL  edge  (ft/s) 
angle  of  attack  (deg) 
altitude  (ft) 

free-stream  speed  of  sound  (ft/s) 

speed  of  sound  at  stagnation  point  (ft/s) 

speed  of  sound  at  Eckert’s  reference  enthalpy  (ft/s) 

speed  of  sound  at  wall  (ft/s) 

velocity  gradient  (ft/s-ft) 

mass  flow  thickness  (ft) 

compressible  skin  friction  coefficient 

incompressible  skin  friction  coefficient 

convergence  identifier 

pressure  coefficient 

specific  heat  of  air  at  constant  pressure  (Btu/lbm-°F) 
mass  flow  thickness  change  (ft) 

BL  thickness  (ft) 

stored  values  of  BL  thickness  (ft) 

displacement  thickness  (ft) 

stored  values  of  BL  displacement  thickness  (ft) 

BL  displacement  thickness  (ft) 
stagnation  point  BL  thickness  (ft) 
momentum  thickness  change  (ft) 

BL  thickness  at  nose/body  tangency  point  (ft) 

incident  angle  at  nosc/body  tangency  point  (deg) 

cone  half  angle  (deg) 

change  in  incidence  angle  (deg) 

body  angle  at  nose/body  tangency  point  (deg) 

geometry  identifier 

momentum  thickness  (ft) 

surface  distance  step  (ft) 

degrees  to  radian  factor  (0.01745) 

free-stream  stagnation  density  (lbm/ft3) 

dummy  variable 

dummy  variable 

velocity  gradient  (ft/s-ft) 

velocity  gradient  at  stagnation  point  (ft/s-ft) 

velocity  change  (ft/s) 

change  in  x  coordinate  of  body  surface  for  von  Karman 
radome  shape  at  angle  of  attack  (ft) 
change  in  y  coordinate  of  body  surface  for  von  Karman 
radome  shape  at  angle  of  attack  (ft) 

BL  edge  Mach  no./free-stream  Mach  no. 

BL  edge  pressure/free-stream  pressure 
sharp  cone-attached  shock  angle  (deg) 
entrainment  rate  (ft/s-ft) 
scale  factor  (currently  1.0) 
free-stream  Mach  no. 
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FSX 

GAM 

GAMMA 

GEO 

H 

HI 

H1K 

HBAR 

HE 

HEX 

HK 

HO 

HR 

HREF 

HRL 

HRT 

HRX 

HS 

HTO 

HTX 

HW 

I 

11-15 

IAIR 

IFLAG 

IFLAG1 

IGE 

IQS 

ISTEP 

ITR 

ITY 

IX 

J 

JSTEP 
JUNK,  1,2,4 
KKK 
LAM 
LT 

MACH 

MARKER 

MDOT 

ME 

ME2 

MMARKER 

MNORM 

MO 

M02 

MRL 

MRX 

MST 

MUE 

MUREF 


RMS*VE*Y2 
ratio  of  specific  heats 
bow  shock  angle  (deg) 
geometry  identifier 

compressible  conventional  shape  parameter  (DD/DM) 
compressible  modified  shape  parameter  (BIGD/DM) 
incompressible  modified  shape  parameter 
transformed  shape  parameter 
BL  edge  enthalpy  based  on  0°R  (Btu/lbm) 

BL  edge  enthalpy  based  on  540°R  (Btu/lbm) 
incompressible  conventional  shape  parameter 
free-stream  static  enthalpy  based  on  0°R  (Btu/lbm) 
recovery  enthalpy  based  on  0°R  (Btu/lbm) 
reference  enthalpy  for  Cohen  correlation  (8465  Btu/lbm) 
last  (previous)  recovery  enthalpy  based  on  O’R  (Btu/lbm) 
recovery  enthalpy  based  on  0*R  (Btu/lbm) 
recovery  enthalpy  based  on  540°R  (Btu/lbm) 

Eckert’s  reference  enmalpy  based  on  0°R  (Btu/lbm) 

total  enthalpy  based  on  0°R  (Btu/lbm) 

total  enthalpy  based  on  540°  R  (Btu/lbm) 

wall  enthalpy  based  on  0°R  (Btu/lbm) 

dummy  indexing  variable 

1-5 

switch  between  perfect  and  equilibrium  air  properties 

no.  of  laminar  streamline  swallowing  iterations  to  converge 

no.  of  times  laminar  streamline  swallowing  iterations  have 

been  hung  up  in  loop  at  a  discontinuity 

geometry  specification  variable 

currently  not  used 

indexes  for  each  flow  station 

transition  criterion  index 

initial  condition  index 

index  for  each  change  in  NPTS 

dummy  indexing  variable 

steps  for  each  flow  case 

dummy  variables 

equals  zero  for  first  iteration 

laminar  BL  identifier 

length  of  radome  to  tangency  point  (ft) 

Mach  no. 

steps  for  each  marching  step 
mass  flow  rate  (lbm/s) 

Mach  no.:  at  BL  edge 

Mach  no.:  at  BL  edge  squared 

equals  1  to  initialize  turbulent  swallowing  equations 

free-stream  Mach  no.  normal  to  bow  shock 

free-stream  Mach  no. 

free-stream  Mach  no.  squared 

last  (previous)  BL  edge  Mach  no. 

interpolated  Mach  no.  at  BL  edge 

viscosity  at  Eckert’s  reference  enthalpy  (lbm/ft-s) 

viscosity  at  BL  edge  (lbm/ft-s) 

reference  viscosity  for  Cohen  correlation  (1.15312E-5  lbm/ft2-s) 
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MUS 

viscosity  at  stagnation  point  (lbm/ft-s) 

NOUT 

output  to  print  and  plot  files  every  NOUT  step 

NPTS 

number  of  stations  wail  temperature  will  be  specified 

NUE 

kinematic  viscosity  at  BL  edge  (ft2/s) 

NURL 

switch  for  creating  output  files 

NUW 

kinematic  viscosity  at  wall  (ft2/s) 

P 

pressure  (lbf/ft2) 

PC 

surface  pressure  (lbf/ft2) 

PE 

pressure  at  BL  edge  (lbf/ft2) 

PEL 

last  (previous)  BL  edge  pressure  (lbf/ft2) 
interpolated  pressure  at  BL  edge  (lbf/ft2) 

PEX 

PI 

3.14159 

PO 

free-stream  static  pressure  (lbf/ft2) 

PR 

Prandtl  no. 

PREF 

conversion  factor  from  atm  to  psf  (2116.224  lbf/ft2) 

PTO 

free-stream  total  pressure  (lbf/ft2) 

PTS 

pressure  at  stagnation  point  (lbf/ft2) 
static  pressure  at  wall  (lbf/fr) 

PW 

Q 

heat  transfer  rate  (Btu/ft2-s) 

QCW 

cold  wall  heat  transfer  rate  based  on  wall  temperature  of 
540°R  (Btu/ft2-s) 

QHW 

hot  wall  heat  transfer  rate  (Btu/ft2-s) 

QSTAGD 

Detra,  Kemp,  Riddell  stagnation  point  heat  transfer  rate 
(Btu/ft2-s) 

QSTAGF 

Fay-Riddell  stagnation  point  heat  transfer  rate  (Btu/ft2-s) 

QSTAGL 

Lees’s  stagnation  point  heat  transfer  rate  (Btu/ft2-s) 

R 

recovery  factor 

RAD 

1  radian  =  57.3 

RADIUS 

radius  of  body  in  body  coordinate  system  (ft) 

RC 

term  in  equation  for  bow  shock  angle 

RES 

Reynolds  no.  based  on  surface  distance 

RET 

Reynolds  no,  based  on  momentum  thickness 

RETING 

incompressible  Reynolds  no. 

RG 

radius  of  ogive  (ft) 

RHOE 

density  at  BL  edge  (lbm/ft3) 

RHOMUE 

density  x  viscosity  at  BL  edge  (lbm2/ft4-s) 
density  x  viscosity  at  wall  (lbm2/ft4-s) 

RHOMUW 

RHOO 

free-stream  density  (lbm/ft3) 

RHOS 

stagnation  point  density  (lbm/ft3) 

RHOST 

density  at  Eckert’s  reference  enthalpy  (lbm/ft3) 

RHOW 

density  at  wall  (lbm/ft3) 

RHREF 

reference  density  for  Cohen  correlation  (0.0803707 
lbm/ft3) 

RM1 

term  in  equation  for  bow  shock  angle 

RME 

density  x  viscosity  at  BL  edge  (lbm2/ft4-s) 

RMER 

dummy  variable 

RMES 

RME  at  stagnation  point  (lbm2/ft4-s) 

RMS 

density  x  viscosity  at  Eckert’s  reference  enthalpy 
(lbm2/ft4-s) 

RMSR 

dummy  variable 

RMSS 

RMS  at  stagnation  point  (lbm2/ft4-s) 

RMWR 

dummy  variable 

RMWS 

density  x  viscosity  at  wall  at  stagnation  point  (lbm2/ft4-s) 
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RN 

nose  radius  (ft) 

RR 

radius  of  radome  base  at  tangency  point  (ft) 

RRR 

y  coordinate  of  radome  base  including  angle  of  attack  (ft) 

RUC 

interpolated  heat  transfer  coefficient  (lbm/ft2-s) 

RUCH1 

heat  transfer  coefficient  (lbm/ft2-s) 

RUCHL 

last  (previous)  heai  transfer  coefficient  (lbm/ft2-s) 

S 

surface  distance  (ft) 

SI 

sharp  cone  distance  to  tangency  point  (ft) 

SEFF 

effective  sharp  cone  surface  distance  (ft) 

SL 

last  (previous)  surface  distance  (ft) 

SMAX 

maximum  surface  distance  (ft) 

SO 

initial  surface  distance  (ft) 

SQRE 

square  root  of  VE  >:•  Sl/NUW 

SR 

y  coordinate  of  body  in  free-stream  coordinate  system  (ft) 

SRTP 

y  coordinate  of  nose/body  tangency  point  in  free-stream 
coordinate  system  (ft) 

ss 

entropy  at  BL  edge  (Btu/lbm-°R) 

sso 

free-stream  entropy  (Btu/lbm-*R) 

SS2 

entropy  at  stagnation  point  (Btu/lbm-*R) 

SST 

entropy  at  Eckert’s  reference  enthalpy  (Btu/lbm-*R) 

SSW 

entropy  at  wall  (Btu/lbm-°R) 

STA 

wall  temperature  (*R) 

STN 

surface  distance  to  nose/body  tangency  point  (ft) 

STRAN 

transition  Reynolds  no. 

TANGAM 

tangent  of  GAMMA 

TDUM 

dummy  variable 

TE 

BL  edge  temperature  (°R) 

TERM  1-9 

terms  in  shape  factor  conversion  equations 

TGJ 

50,063  (lbm-ft2/Btu-s2) 

TH 

body  incidence  angle  to  flow  (deg) 

TH1 

body  slope  at  nose/body  tangency  point  with  respect  to 
free-stream  (deg) 

THL 

last  (previous)  incidence  angle  (deg) 

T1 

time  (s) 

TIME 

time  (s) 

TO 

free-stream  static  temperature  (°R) 

TOO 

total  temperature  (°R) 

TR 

recovery  temperature  (°R) 

TRAN 

specifies  transition  correlation 

TST 

temperature  at  Eckert’s  reference  enthalpy  (°R) 

TTO 

free-stream  total  temperature  (°R) 

ITS 

stagnation  point  temperature  (°R) 

TURB 

turbulent  BL  identifier 

TYPE 

specifies  format  of  input  conditions 

V2 

velocity  immediately  behind  shock  (ft/s) 

VE 

velocity  at  BL  edge  (ft/s) 

VE1.2 

saved  BL  edge  velocities  (ft/s) 

VENUW 

VE  X  NUW  (ft  Vs2) 

VNORM 

free-stream  velocity  normal  to  bow  shock  (ft/s) 

VO 

free-stream  velocity  (ft/s) 

VTAN 

fret  stream  velocity  tangential  to  bow  shock  (ft/s) 

WE 

molecular  weight  at  BL  edge  (lbm/ibm-mol) 

WO 

free-stream  molecular  weight  (lbm/lbm  mol) 
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ws 

WST 

WW 

X 

XLEN 

XT 

XTC 

XXXI, 2 

Y 

YS 

YT 

YTC 

YYY1,2 

ZO 


molecular  weight  at  stagnation  point  (lbm/lbm-mol) 
molecular  weight  at  Eckert’s  reference  enthalpy 
(lbm/lbm-mol) 

molecular  weight  of  gas  at  wall  (lbm/lbm-mol) 
x  coordinate  of  body  in  free-stream  coordinate  system  (ft) 
surface  distance  (ft) 

x  coordinate  of  nose/body  tangency  point  at  angle  of  at¬ 
tack  in  free-stream  coordinate  system  (ft) 
surface  x  coordinate  in  body  coordinate  system  at  nose/ 
body  tangency  point  (ft) 

x  coordinate  of  body  surface  for  von  Karman  radome  shape 
at  angle  of  attack  (ft) 

y  coordinate  of  body  in  free-stream  coordinate  system  (ft) 
radius  of  free-stream  stream  tube  required  to  match  MDOT 
(ft) 

y  coordinate  of  nose/body  tangency  point  at  angle  of  at¬ 
tack  in  free-stream  coordinate  system  (ft) 
surface  radius  in  body  coordinate  system  at  nose/body  tan¬ 
gency  point  (ft) 

y  coordinate  of  body  surface  for  von  Karman  radome  shape 
at  angle  of  attack  (ft) 
altitude  (ft) 
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GLOSSARY 


Cf 

-  skin  friction 

T 

-  temperature 

F 

-  entrainment  parameter 

TT 

-  total  temperature 

h 

-  enthalpy 

TW 

-  wall  temperature 

H 

-  conventional  form 

u 

-  velocity 

parameter  (5*/0) 

X 

-  longitudinal  coordinate 

Hi 

-  modified  form  parameter 

y 

-  radius  from  centerline 

(S-SV0) 

a 

-  angle  of  attack 

m 

-  boundary  layer  mass  flow 

0 

-  velocity  gradient  ( dut  /ds) 

rate 

7 

-  ratio  of  specific  heats  ( cp/cv ) 

M 

-  Mach  no. 

6 

-  boundary  layer  thickness 

P 

-  pressure 

6* 

-  boundary  layer  displacement 

Pr 

-  Prandtl  no. 

thickness 

PT 

-  total  pressure 

A 

-  mass  flow  thickness 

Q 

-  heat  transfer  rate 

V 

-  sharp  cone-attached  shock 

QS 

-  stagnation  point  heat  trans¬ 

angle 

fer  rate 

6 

-  boundary  layer  momentum 

r 

-  radius  of  body  at  S 

thickness 

R 

-  radius 

A 

-  bow  shock  angle 

Re 

-  Reynolds  no. 

-  viscosity 

S 

-  surface  distance  from  stag¬ 

V 

-  kinematic  viscosity 

nation  point 

P 

-  density 

Sy 

-  function  defined  in  Ref.  26 

a 

-  cone  incident  angle  to  flow 

Subscripts 

c  -  sharp  cone 

sh 

-  shock 

e  -  boundary  layer  edge 

w 

-  wall 

n  -  blunted  nose 

* 

-  Eckert’s  reference  temperature 

r  -  recovery 

s  -  stagnation  point 

00 

-  free  stream 
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